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1.  INTRODUCTION 


"A  workshop  on  Natural  Charging  of  Large  Space  Structures  in 
Near  Earth  Polar  Orbits,  sponsored  by  the  Air  Force  Geophysics 
Laboratory  and  Boston  College,  was  held  on  14  and  15 
September,  1982,  at  AFGL.  Many  of  the  specialists  in  the  area  of 
spacecraft  charging  problems  assembled  to  discuss  the  necessity 
and/or  possibility  of  developing  realistic,  effective  codes 
describing  the  interaction  of  low-earth  orbiting  systems  with  the 
environment.  ...There  are  considerable  differences  of  opinion 
among  specialists  in  this  area  concerning  the  nature  and 
magnitude  of  the  problems  and  the  feasibility  of  developing 
realistic  models  or  codes  to  deal  with  them. " 1 

The  POLAR  code  was  started  amidst  controversy.  At  the  workshop,  questions 
were  raised  about  whether  large  objects  would  charge  in  polar  orbit  and 
whether  it  was  possible  to  model  the  relevant  physics.  Others  expressed 
concerns  about  how  resultant  computer  models  could  be  validated.  The  polar 
orbit  charging  question  was  answered  and  observed  by  Geophysics  Laboratory 
(GL)  instruments  mounted  on  DMSP  spacecraft. 


The  POLAR  code  was  written  to  model  charging  by  large  spacecraft  in  low, 
polar  orbit.  This  report  documents  comparisons  of  POLAR  code  calculations 
with  flight  experiments.  Calculations  of  the  plasma  wake  behind  the  Shuttle 
Orbiter  are  compared  with  in  the  situ  measurements  of  Murphy  et  al.  [1989]. 
Calculations  of  the  charging  of  DMSP-7,  performed  by  Dr.  David  Cooke  of  the 
GL,  are  compared  with  the  flight  measurements  of  Gussenhoven  et  al.  [1985]. 
Calculations  of  current  and  potential  on  the  SPEAR  I  rocket  are  compared  with 
the  flight  measurements.  Calculation  of  electron  collection  by  the  CHARGE  II 
mother  rocket  payload  are  compared  with  the  flight  observations  of  Myers  et  al. 
[1989].  In  all  cases,  the  POLAR  calculations  agreed  with  the  essential  features 
of  the  observations. 


The  Shuttle  Orbiter  and  SPEAR  I  comparisons  were  intended  to  be  with 
preflight  predictions.  In  both  cases,  however,  malfunctions  during  the  flights 

1  Proceedings  of  the  Air  Force  Geophysics  Laboratory  Workshop  on  Natural  Charging  ot  Large 
Space  Structures  in  Near  Earth  Polar  Orbit  .  14-15  September  1982. 
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changed  important  parameters  and  required  new  calculations.  The  Orbiter 
main  engine  shut  down  early,  lowering  the  orbit  by  a  hundred  kilometers.  A  cap 
on  the  SPEAR  I  plasma  contactor  didnt  open  allowing  potentials  of  thousands 
of  volts  to  appear  on  the  rocket  chassis.  In  both  cases,  the  POLAR  calculations 
were  redone  using  the  same  model  used  for  the  preflight  calculations,  but  with 
parameters  appropriate  for  the  actual  flight. 

This  validation  document  has  three  parts  in  addition  to  the  introduction:  the  first 
is  a  summary  of  the  physics  learned  during  the  POLAR  development  and 
validation;  the  second  is  a  discussion  of  the  physical  processes  modeled  by  the 
POLAR  code;  the  third  part  is  a  collection  of  papers  that  discuss  the 
comparisons  with  the  flight  data. 
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2. 


SUMMARY  OF  POLAR  ACHIEVEMENTS 


"In  brief,  we  require  of  new  theories  that  they  account  for  the  results 
predicted  by  existing  theories  and  that,  in  addition,  they  predict  new 
results. " 2 

The  POLAR  development  and  validation  efforts  have  lead  to  predictions  and  a 
new  level  of  understanding  about  the  physical  mechanisms  that  determine 
spacecraft  plasma  interactions  with  the  ionosphere.  POLAR  has  answered 
many  long,  outstanding  questions.  The  details  are  contained  in  the  attached 
papers,  and  the  validity  of  the  results  is  born  out  by  direct  comparison  with 
observations.  The  following  paragraphs  summarize  POLAR's  most  significant 
contributions  to  spacecraft  interaction  physics.  In  each  case,  the  results  would 
not  have  been  possible  without  the  data  and  cooperation  of  the  flight 
experimenters. 

i .  The  physical  processes  controlling  the  plasma  wake  of  large  (~10  meter) 
spacecraft  have  been  elucidated  and  compared  with  observation. 

Ions  are  accelerated  into  the  plasma  void  behind  large  satellites  by  density 
gradient  electric  fields.  The  POLAR  comparison  with  in  situ  measurements 
of  the  orbiter  wake  [Murphy  &  Katz,  1989]  showed  that  the  density  over 
shoots  and  wake  structure  seen  in  laboratory  experiments  were  not 
relevant  to  spacecraft  in  the  ionosphere. 

V  POLAR  includes  the  physics  necessary  to  model  wakes. 

2  Auroral  charging  is  predictable  and  increases  with  spacecraft  size. 

POLAR  was  developed  to  pursue  the  speculations  of  Parks  &  Katz  [1981]. 
POLAR  was  just  one  part  of  a  Geophysics  Laboratory  (GL)  investigation 
into  auroral  charging.  Measurements  of  auroral  charging  on  DMSP  were 
reported  by  Gussenhoven  et  al.  [1985]  on  DMSP.  POLAR  calculations  by 

2  M  S.  Gussenhoven,  Proceedings  of  the  Air  Force  Geophysics  Laboratory  Workshop  on 

Natural  Charging  of  Large  Space  Structures  in  Near  Earth  Polar  Orbit;  14-15  September  1982, 
Chapter  31 ,  "Requirements  for  Validating  System  Models". 
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Cooke  et  al.  [1989]  using  data  provided  by  Gussenhoven,  reproduced  the 
observed  levels  of  charging.  The  POLAR  calculations  also  showed  that  a 
small  satellite  would  have  charged  to  just  a  few  volts. 

V  POLAR  includes  the  physics  necessary  to  model  auroral  charging. 

3.  Electron  collection  from  the  ionosphere  by  active  experiments  is 
magnetically  limited  except  at  low  (<250  km)  altitudes. 

POLAR  calculations  of  the  CHARGE  II  mother  payload  showed  that  the 
data  of  Meyers  et  al.  [1989]  actually  agreed  with  magnetic  limiting,  contrary 
to  much  speculation  in  the  literature.  The  POLAR  model  of  the  rocket  had 
the  correct  finite  cylinder  geometry,  and  collected  current  in  agreement 
with  the  measurements.  The  POLAR  currents  were  substantially  larger 
than  the  spherical  probe  results  previously  applied  to  the  problem.  The 
SPEAR  I  calculations  [Katz  et.  al.,  1989]  showed  how  deviations  from 
spherical  symmetry  enable  electrons  to  cross  magnetic  field  lines  and  be 
collected  by  high-voltage  probes.  In  neither  case  was  any  "turbulence" 
necessary  to  explain  the  experimental  results. 

V  POLAR  includes  the  physics  necessary  to  model  plasma  sheaths. 
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3. 


POLAR  CODE  PHYSICAL  MODELS 


Below  are  outlined  the  basic  physics  models  solved  in  POLAR.  More  details  on 
the  algorithms  used  to  implement  the  models  are  contained  in  the  POLAR 
Code  Users  Manuals. 

3.1  PLASMA  SHEATH 

POLAR  solves  Poisson's  equation  in  the  volume  surrounding  a  spacecraft. 

2  P 

V0  =  f-  (i) 

o 

For  potentials,  p,  small  compared  with  the  background  plasma  temperature  0, 
the  space-charge  density,  p  ,  is  calculated  assuming  linear  screening: 

~  =  '1^1  ~  ^  sheath  &) 

where  pi-hQ..u  =  0ln2,  (3) 

the  potential  where  linear  shielding  breaks  down  [Parrot,  Storey,  Laframboise  & 
Parker,  1982].  In  many  cases,  the  value  used  in  POLAR  for  Psheath  is  higher, 
due  to  numerical  accuracy  considerations,  but  even  for  the  most  extreme  cases, 
it  is  less  than  a  few  volts.  This  value  sets  the  the  low  potential  accuracy  limit  of 
the  code. 

For  larger  potentials,  the  space  charge  density  is  calculated  by  tracking 
representative  macro  particles  in  a  fixed  potential, 

x  =  ^(x  x  B  -  Vp).  (4) 

The  potentials  are  assumed  to  remain  steady  during  a  particle  transit  through 
the  sheath.  This  sets  the  lower  bound  on  the  timescales,  which  can  be  modeled 
using  POLAR.  Potentials  and  charge  densities  are  iterated  until  a  self- 
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consistent  solution  is  obtained.  To  speed  convergence,  an  analytical  estimate 
of  the  space-charge  density  is  used  to  obtain  the  initial  potential  iterate. 

3.2  NATURAL  CHARGING 

Surface  potential,  <j>,  changes  are  calculated  using  the  capacitance,  C,  between 
the  surface  and  the  underlying  chassis  ground, 

dp 

C‘dT=jne,  XA-  <5) 

where  A  is  the  surface  area.  The  net  current  density  has  several  components. 

jnet  -  jenCident  X  (l  -  Yt*COndar*  a  -  B^ackscat,ef  e)  +  j&aa,h  X  (l  +  Yf^"^  e)  -  jPhoto  (6) 

The  incident  ion  currents  are  calculated  by  tracking  particles  through  the 
sheath.  The  incident  electron  current  has  both  a  low  and  high  energy 
component.  The  low  energy  component  is  treated  as  a  single  temperature 
Maxwellian  and  accurately  represents  the  cold,  dense  Ionospheric  background. 
The  high  energy  component  is  used  to  represent  auroral  electron  fluxes  and 
can  be  of  the  form 

l(E)  =  aE(E-e$)''“*”  +  nVi;  (-J-r')Ee^,,l'*'+bEe  "  *  (7). 

V*<kT)  ) 

This  form  has  been  shown  to  fit  inverted-V  spectra  [Fontheim  et  al.,1982]. 

3.3  PLASMA  WAKE 

The  model  of  the  wake  structure  used  by  POLAR  depends  on  the  position 
relative  to  the  so-called  ion  front.  This  ion  front  marks  the  boundary  where 
electron  density  begins  to  change  on  a  scale  commensurate  with  the  Debye 
length  and  the  ion  density  takes  a  sudden  and  dramatic  drop.  Ahead  of  the  ion 
front,  the  plasma  is  treated  as  rarefied;  its  motion  is  controlled  by  the  thermal 
spread  in  ion  velocities.  Behind  the  front,  the  motion  is  controlled  by  the 
electron  temperature  and  ion  mass. 
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The  governing  equations  in  this  region  ahead  of  the  front,  considering  that 
electrons  are  more  mobile  than  ions  and  that  they  maintain  equilibrium  with  a 
local  potential,  are: 


the  Boltzman  relation; 

ne  =  n0exp(e  d/kTe) 

(8) 

3n,  0 

continuity; 

nr +  "37  tn'v) " 0 

(9) 

equation  of  motion; 

f)v  vr)v  -  e 
<A  +  f)z  ~  M  Dz 

(10) 

Poisson's  equation; 

dz2  =4rte(ne-n) 

(11) 

where  n0 

ambient  density 

n, 

ion  density 

electron  density 

ee  = 

electron  temperature 

e 

electron  charge 

0 

local  potential 

k 

Boltzman's  constant 

z  is  a  variable  representing  distance  parallel  to  the  front  velocity  or,  in  the  case 
of  POLAR,  perpendicular  to  the  orbital  velocity. 

Crow  et  al.  [1975]  have  numerically  solved  the  above  equations  to  predict  the 
position  of  the  ion  front.  POLAR  uses  an  analytical  fit  to  the  Crow  results  [Katz  et 
al.,  1985]: 


p(t)  =  2X  d  {  (  a>t  +  In  ( 1  +  aoit)  -  tot  -  ( 1-^jp)  (cot  -  ^  ln(  1  +  a  cot))  j* 


where  to  = 


(47tnoe?)  (k09) 


M 


= 


47rn0e 


(12) 
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are  the  ion  plasma  frequency  and  Debye  length,  a  is  a  free  parameter 
determined  to  be  ~  1 .6. 


Katz  et  al.  [1985]  showed  that  this  formula  agrees  well  with  laboratory  data  from 
Wright  et  al.  [1985].  Ahead  of  this  front  zF,  the  plasma  is  assumed  to  expand 
due  to  thermal  motion,  the  so-called  "neutral  approximation".  Behind  zF  the 
plasma  evolves  into  a  state  that  is  self-similar  [Chan  et  al.,  1984],  The  self¬ 
similar  solution  of  equations  8-11  for  z  >  S0t  is 


n-n-“pH£s71)! 


where  S0 


(kee/M)2 


is  the  ion  acoustic  speed. 


We  take  the  time  variable  to  be  spatially  defined  as: 


(13) 


(14) 


where  x  is  the  distance  behind  the  object  (perpendicular  to  z)  and  V0  is  the 
orbital  velocity. 

The  wake  routines  in  POLAR  employ  both  limiting  cases:  (1)  ahead  of  the  ion 
front  the  electric  field  is  negligible  and  the  motion  of  ions  is  identical  to  neutrals; 
(2)  behind  the  ion  front,  whose  position  is  determined  by  equation  (12),  the 
quasi-neutral,  self-similar  solution  of  equation  (13)  is  implemented. 


POLAR  has  routines  that  accurately  model  the  geometry  of  the  object,  and  the 
"neutral  ion"  trajectories  are  calculated  from: 


fjQU)  =  g(*.n)  fto  (?)  (15) 

where  f  (tf)is  the  unperturbed  distribution  function  for  a  drifting  Maxwellian; 
g(2,fi)  has  value  "0"  if  a  ray  starting  from  9  and  going  in  the  direction  LI  would 
strike  the  vehicle,  and  "1"  if  it  would  not. 


8 


The  local  density  is  given  by: 


n,  =  (*)J  r  (*.v)  =  Jg(*,n)  {jfto(yfQ)v2dv}dn  (16) 

This  initial  density  calculated  in  three  dimensions  for  neutral  particles  is 
compared  with  density  calculated  assuming  the  complex  geometric  object  is 
replaced  by  a  flat  plate  at  a  position  where  the  dominant  source  appears  at  the 
object  edge.  This  ratio  provides  a  "geometric  correction  factor"  that  is  applied  to 
the  quasi-neutral,  one-dimensional  solution  discussed  above  for  positions 
behind  zF.  In  this  way,  POLAR  can  calculate  quite  rapidly  an  approximate  value 
for  the  ion  and  electron  densities  in  the  wakes  of  complex  objects. 

Note  that  the  assumptions  behind  the  front  are:  (1)  the  electron  temperature 
and  ion  mass  govern  the  equation  of  motion;  (2)  the  plasma  is  quasi-neutral; 

(3)  the  magnetic  field  does  not  affect  the  ion  or  electron  motion;  (4)  equa¬ 
tion  (12)  serves  as  a  good  approximation  for  determining  the  boundary  of  the 
ion  front;  and  (5)  the  geometric  correction  factor  calculated  in  detail  with  the 
3-dimensional  neutral  model  can  be  approximately  applied  to  correct  the 
plasma  densities  as  well.  Therefore,  the  algorithm  can  address  complex 
geometries  but  takes  advantage  of  the  smooth  wake  structure  characteristic  of 
ionospheric  plasmas  where  0j/0e  =  1. 
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4.  FLIGHT  DATA  COMPARISON  PAPERS 


The  papers  on  the  following  pages  discuss  POLAR  code  model  comparisons 
with  flight  data. 


JOI  RNAI  Of-  GFOPHS  SIC'Al.  RFSEARCH.  VOL  V4.  NO  A2.  PAGES  1450-I4SK,  FEBRUARY  |.|4H9 


Structure  of  the  Bipolar  Plasma  Sheath  Generated  by  SPEAR  I 

I  kxr/.'  G  A  Josgfwarp.1  V  A  Davis,1  M  J.  Mandfi.i.1  R  A,  Kuharski.1  J.  R,  Lilley.  Jr.,1  W.  J.  Raitt,*’ 

D  L  Cooke.'  R  B  Torbfrt.4  G  Larson.4  and  D  Rad4 

I  he  Spate  Power  Experiment  Aboard  Rockets  I  (SPEAR  I)  biased  two  10-cm  radius  spheres  as  high 
js  46.UOO  V  positive  with  respect  to  an  aluminum  rocket  body  The  experiment  measured  the  steady  state 
current  to  the  spheres  and  the  float, ng  potential  of  the  rocket  body  Three-dimensional  calculations 
performed  using  N  ASC  AP  LEO  and  POLAR  2  0  show  that  both  ion-collecting  and  electron-collecting 
sheaths  were  formed  The  rocket  body  potential  with  respect  to  the  ionospheric  plasma  adjusted  to 
achieve  a  balance  between  the  electron  current  collected  by  the  spheres  and  the  secondary  electron- 
enhanced  ion  current  to  the  rocket  body  This  current  balance  was  obtained  with  a  large  ton-collecting 
sheath  that  enveloped  most  of  the  electron-collecting  sheath  and  reduced  the  area  for  collection  of 
ionospheric  electrons  The  calculated  current  is  in  agreement  with  the  flight  measurement  of  a  steady 
State  current  of  less  than  I  tl)  A  The  calculations  show  that  the  rocket  body  was  driven  thousands  of 
volts  negative  with  respect  to  the  ionospheric  plasma  The  calculated  rocket  potential  is  within  the 
uncertainly  of  that  inferred  from  ion  spectrometer  data  The  current  flowed  through  the  space  plasma 
There  was  almost  no  direct  charge  transport  between  the  spheres  and  the  rocket  body. 


iMRODLt  won  Power  Experiment  Aboard  Rockets  I  (SPEAR  I)  (W  J  Raitt 

Since  the  first  time  that  electron  guns  were  placed  on  rock-  al .  unpublished  manuscript.  19X8)  was  specifically  designed 
ets.  the  voltage  on  a  rocket  necessary  to  collect  ionospheric  10  measure  whether  or  not  the  Earth's  magnetic  field  impedes 
electrons  and  complete  the  circuit  has  been  the  subject  of  electron  collection  SPEAR  I  had  two  10-cm  radius  spherical 
much  debate  [MmrUer.  1980]  Classical  theories  f Langmuir  electron  collectors  each  at  the  end  of  a  3-m  boom  The 
11,1,1  Hlodqett.  1924:  Unit-Smith  and  l.anqmutr.  1926:  Beard  planned  experiment  w3s  to  apply  a  positive  voltage  on  the 
.1111/  Jnhnnnt.  1961.  (  lien.  1965.  l.afrunibotse  and  Rubinstein,  spheres,  up  to  46.000  V  with  respect  to  the  rocket  body,  and 
I9~6.  Rubinstein  and  l.alrambmse  1978.  1982.  1983.  Parker  measure  the  current  of  collected  electrons  from  the  lono- 
,md  Uurphx.  1967]  predict  that  even  for  low-current  electron  sphere  In  place  of  an  electron  beam.  SPEAR  I  carried  a 
beams  (less  than  100  mAl.  rocket  bodies  will  charge  to  thou-  plasma  contactor  on  the  bottom  of  the  rocket  to  emit  elec- 
sands  of  volts  Early  experiments,  however,  measured  poten-  |fons  and  complete  the  circuit  through  the  ionosphere.  The 
Hals  of  less  than  100  V  for  beam  currents  as  high  as  an  ampere  Plasma  contactor  was  designed  to  keep  the  rocket  body  close 
tsee  Winkler  [1980]  for  a  review  of  the  experiments).  Theor-  IO  ,*ie  Siimc  potential  as  the  surrounding  ionosphere.  Conse- 
etictans  questioned  the  validity  of  the  measurements;  expert-  *-l uen 1 1>.  most  of  the  potential  difference  applied  between  the 
ments  dismissed  the  theories  as  obviously  irrelevant.  In  more  *Phere  and  rocket  body  was  intended  to  be  dropped  across  a 
recent  years,  experiments  performed  above  300  km  produced  large  electron-collecting  sheath. 

measurements  more  in  agreement  with  theoretical  predictions  During  the  SPEAR  I  flight,  the  cover  failed  to  come  off  the 
[  \laehtum  el  al  .  1988]  Using  tethered  subsatelliles  (so-called  P'asma  contactor.  Because  of  this,  the  plasma  contactor  was 
mother  daughter  "  payloads),  potentials  have  been  measured  uni|bie  to  expel  the  electrons  collected  by  the  spheres.  The 
more  reliably  While  at  low  altitudes  the  measured  rocket  po-  r°cket  body  was  driven  thousands  of  volts  negative  to  achieve 
tenttals  remain  low.  the  results  for  higher-altitude  rockets  lend  current  balance  Instead  of  a  single,  nearly  spherical  electron- 
to  agree  with  classical  theories  of  space-charge  limited  col-  collecting  sheath  around  the  spheres,  a  second,  larger  ion- 

lectton  The  low -altitude  results  are  thought  to  stem  from  ioni-  collecting  sheath  was  formed  around  the  rocket  body,  par- 

ration  of  neutral  gas  within  the  electron  collecting  sheath  [Lai  Hally  enveloping  the  electron  sheath.  The  structures  and  inter¬ 
im  al .  1985]  Still  unresolved,  ow  ing  to  the  relatively  low  po-  actions  between  the  two  sheaths  and  the  closure  of  the  circuit 
tenttals  on  the  rockets,  ts  the  influence  of  the  Earth's  magnetic  through  the  ionosphere  are  the  subjects  of  this  paper.  The 
held  on  electron  collection  basic  theory  of  high-voltage  objects  immersed  in  a  plasma  is 

L  lassical  theories  predict  that  near  large  high-voltage  presented,  followed  by  the  results  from  detailed,  three- 
probes.  the  ionospheric  electrons  are  constrained  more  by  dimensional  calculations  for  the  SPEAR  I  parameters.  Both 
magnetic  field  lines  than  attracted  by  the  probe  potential.  This  *lw  flight  results  and  the  supporting  calculations  show  that 

should  result  tn  a  decrease  m  current  from  space -charge  limit-  when  the  neutral  gas  pressure  is  low  enough  that  ionization 

ed  predictions  that  has  yet  to  be  observed  in  space  The  Space  van  be  ignored,  the  sheath  currents  are  small  and  are  con¬ 
trolled  by  the  geometry  of  the  exposed  high-voltage  compo¬ 
nents 


EXHtRIMFNTAt  ARRANGFMFNT 

The  geometric  arrangement  of  the  booms  and  payload 
structure  is  shown  in  Figure  I.  The  10-cm  radius  gold-plated 
spheres  were  mounted  on  bushings  constructed  with  grading 
rings  connected  by  resistors  The  bushings  produced  a  uniform 
potential  gradient  from  the  spheres  to  payload  ground  at  the 
Y  junction  of  the  two  bushings  with  the  plastic  support  boom 
shown  in  figure  I  The  lotal  resistance  along  each  bushing 
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Surtac*  Material*  Dfmentioni  (matart) 


fig  I  The  SPLAR  I  experimenl  The  rocket  body  was  alu¬ 
minum.  ihe  booms  were  plastic  the  bushings  were  nickel,  and  the 
spheres  were  gold 

was  I  Mfl,  and  (he  number  of  rings  used  resulted  in  a  poten¬ 
tial  drop  of  about  I  kV  between  each  pair  of  rings  when  the 
sphere  was  at  a  potential  of  46  kV.  The  rocket  body  was 
polished  aluminum  with  a  surface  area  of  just  over  5  nr  and 
was  connected  to  payload  ground. 

The  voltage  between  the  spheres  and  the  payload  ground 
was  supplied  from  a  charged  2.5-/rF  capacitor  switched  to  the 
spheres  by  a  high-voltage  relay.  A  current-limiting  resistor  of  2 
kO  was  connected  in  series  with  the  capacitor.  The  electric 
schematic  diagram  of  the  high-voltage  section  for  one  sphere 
is  shown  in  Figure  2.  The  fixed  resistor  of  700  kQ  in  parallel 
with  the  bushing  resistors  of  1  Mfl  resulted  in  a  decay  time 
constant  of  I  s  for  the  capacitor,  ensuring  that  even  in  the 


absence  of  plasma  currents,  the  applied  potential  would  decay 
to  near  0  V  during  the  5-s  period  for  which  the  capacitor  was 
connected  lo  the  sphere. 

Figure  2  illustrates  the  two  current  monitors  used  to  mea¬ 
sure  the  total  current  from  the  capacitor  and  the  current  flow¬ 
ing  down  the  graded  bushing  The  plasma  current  to  the 
spheres  was  determined  from  the  relation 


where  F  is  the  measured  voltage  on  the  capacitor  in  kilovolts 
and  the  constant  700  is  the  constant  fixed  resistance  connected 
between  the  capacitor  and  payload  ground  An  example  of  the 
plasma  current  to  the  sphere  as  a  function  of  applied  voltage 
is  shown  in  Figure  .V 

In  addition  to  the  high-voltage  generating  and  monitoring 
instruments,  the  payload  also  contained  diagnostics  to  moni¬ 
tor  neuttal  pressure,  ambient  electron  density  and  temper¬ 
ature.  VLF  and  HF  wave  intensities,  energetic  ion  and  elec¬ 
tron  fluxes  up  to  JO  keV.  and  optical  imagery  and  photometry 
of  glow  in  the  vicinity  of  the  spheres.  A  hollow  cathode 
plasma  contactor  was  also  included  in  the  payload  with  its 
plasma  orifice  at  the  end  of  the  payload  removed  from  the  HV 
spheres. 

SPFAR  I  was  launched  from  NASA  Wallops  Flight  Facility 
on  December  I  T.  19X7.  at  2045  EST  It  reached  an  apogee  of 
269  km  at  251  s  after  launch.  All  of  the  results  shown  in  this 
paper  were  taken  within  10  km  of  350  km  No  significant 
altitude  effect  was  observed  between  120  and  320  km.  The 
pressure  measured  during  most  of  the  high-voltage  experi¬ 
ments  ranged  between  10  *  and  10' 5  tort  (1  Pickett  and  G 
Murphy,  private  communication.  19881  This  is  significantly 
higher  than  the  unperturbed  ambient  environment.  During  the 
last  two  high-voltage  shots.  Ihe  rocket  was  below  130  km  and 
the  background  pressure  was  in  excess  of  10  *  torr.  During 
these  last  shots,  extensive  iom/atton  occurred  and  our  calcula- 
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Potential  (kV) 

f  ; jz  i  The  current  versus  volume  Jala  where  only  one  sphere  was 
charged  The^e  data  were  typical  of  all  the  higher-altitude  shols. 

nonc  Jo  not  apply  For  most  of  the  experiment,  the  rocket  was 
oriented  so  that  the  Earth's  magnetic  held  was  perpendicular 
to  the  plane  of  the  bushings  and  booms.  During  the  later 
portion  of  the  flight,  the  rocket  was  reoriented  with  respect  to 
the  magnetic  field,  but  no  effect  was  seen  :n  the  data. 

The  failure  of  (he  cap  remoyal  on  the  plasma  contactor 
resulted  in  the  payload  ground  being  driven  to  high  negative 
potentials  as  described  earlier  The  resulting  ion  flux  from  the 
ionospheric  plasma  to  the  vehicle  was  detected  by  the  ion 
sensors  of  the  energetic  particle  detectors.  The  ion  flux  energy 


spectra  showed  a  high-energy  cutoff  which  varied  as  the  ap 
plied  vehicle-sphere  potential  varied.  The  expected  peak  in  the 
ion  spectra  was  observed,  but  it  was  broader  than  expected 
from  classical.  colliMonless  acceleration  of  ionospheric  tors 
through  ihe  vehicle  charge  sheath.  Upper  and  lower  bounds 
on  the  vehicle  potential  based  on  the  spectral  cutoff  and  peak 
energies  are  shown  in  Figure  4. 

F  urther  details  on  the  instrumentation  and  examples  of  the 
flight  data  achieved  can  be  found  elsewhere  (W  J  Raitt  et  al . 
unpublished  manuscript.  I9XKI 

Thiory 

1  he  calculations  presented  depend  on  a  few  basic  principles 
First,  a  space  charge  sheath  forms  around  the  rocket  surfaces 
which  are  at  high  potentials  with  respect  to  the  surrounding 
ionospheric  plasma  The  size  of  the  sheath  i»  such  that  the  net 
space  charge  in  the  sheath  cancels  the  charge  on  the  surfaces. 
Second,  rocket  components  positive  with  respect  to  the  iono¬ 
spheric  plasma  collect  electrons,  and  negative  surfaces  both 
collect  ions  and  emit  secondary  electrons.  Finally,  charges 
flow  between  the  SPF.AR  I  rocket  and  the  ionosphere  in  a 
complete  circuit :  the  sum  of  the  current  from  the  positive 
surfaces  is  exactly  balanced  by  the  current  to  the  negative 
surfaces  These  basic  principles  control  the  lime-averaged  re¬ 
sponse  between  SPF.AR  I  and  the  ionosphere 

A  high-voltage  capacitor  in  the  SPEAR  I  rocket  applied  as 
much  as  46.000  V  between  a  10-cm  radius  sphere  and  the 
body  of  the  rocket  In  a  vacuum,  applying  this  potential  to  an 
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t  ig  4  I  Ik*  cspermvemallv  measured  spacecraft  ground  potential  as  .1  function  of  lime  during  one  discharge  of  the 
capacitors  VXSC'AP  l.FO-predicled  ground  potentials  are  shown  jc  diamonds  on  the  curve 
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isolated  sphere  would  cause  a  charge  of  5  x  10  C  lo  appear 
on  (he  surface  If  immersed  in  a  n  =  5  *  10"’  m  ’.  k  lr  =  0  I 
eV  plasma,  typical  of  (he  ionosphere  observed  during  the 
SPEAR  I  llight.  the  potential  excludes  ions  from  a  substantial 
volume  around  the  sphere  If  the  electrons  were  immobile, 
ions  would  be  excluded  from  a  2.5-m  radius  volume  sur¬ 
rounding  the  sphere  ((4ti  3)R'»ii.e  -  4m:(lr<pl  The  space  charge 
of  the  remaining  electrons  would  balance  the  surface  charge 
on  the  sphere  The  transition  from  a  quasi-neutral  plasma  to 
the  complete  exclusion  of  ions  occurs  in  a  few  Debye  lengths 
at  the  edge  of  this  volume.  On  the  scale  of  the  volume  radius, 
this  is  a  very  rapid  transition  f  or  the  ion  space  charge  to  be 
significant,  the  local  potential  must  be  within  a  few  A  T.  of  the 
ambient  plasma  potential.  This  is  less  than  10  4  of  the  applied 
potential  during  the  SPEAR  I  flight  I'nder  these  conditions 
the  sharp  sheath  edge  approximation,  i.e..  the  assumption  that 
the  transition  between  the  quasi-neutral  ambient  plasma  and 
the  ion-free  electron  sheath  is  discontinuous,  is  accurate.  This 
approximation  is  used  in  all  of  the  calculations  below 

When  the  electron  response  is  included,  the  sheath  radius 
expands.  The  self-consistent,  space-charge-limited  electron 
density  and  potential  in  the  sheath  surrounding  a  high-voltage 
spherical  probe  was  first  reported  by  l.tini/nnnr  and  Blodt/ett 
[1924],  The  potential  is  obtained  by  solving  Poisson's  equa¬ 
tion  with  the  boundary  conditions 

01  u  I  =  I 

0|R,)-<)  (2) 

1  r  lr  «. 

where  a  is  the  sphere  radius  and  R%.  the  sheath  radius,  is  found 
as  an  eigenvalue.  The  last  equation  results  from  Gauss's  law 
w  hen  the  space  charge  cancels  the  surface  charge 

Plasma  electrons  enter  at  the  sheath  edge  and  are  absorbed 
by  the  sphere.  In  steady  state  and  without  a  magnetic  field, 
current  continuity  and  conservation  of  energy  lead  to  an  ana¬ 
lytic  expression  for  the  space  charge  I  his  expression  depends 
only  on  the  radial  position,  the  potential  at  that  position,  the 
sheath  radius,  and  the  current  density  entering  she  sheath. 


Parrot  el  til  [  1982]  report  that  weak,  long-range  potentials  in 
the  plasma  increase  the  current  density  entering  the  sheath  to 

.  HI 

where  T,h,rmj|  is  ihc  one-sided  plasma  electron  current  density 
This  sheath  current  results  when  magnetic  fields  are  neglected 
No  comparable  results  are  available  with  magnetic  fields  The 
calculations  use  |4)  to  relate  the  sheath  edge  current  density  to 
the  nlasma  thermal  current.  The  solution  to  Poisson's  equa¬ 
tion  with  the  boundary  conditions  in  (2)  and  the  space  charge 
density  in  (3)  gives  a  sheath  radius  of  7. 1  m 

The  identical  charge  density  could  be  found  by  solving  the 
equations  of  motion  for  the  electrons  and  integrating  over  the 
local  distribution  function  The  analytical  expression  for  space 
charge  density  [Mandell  and  Katz.  1983]  is  simpler  to  evalu¬ 
ate  than  integrating  the  equations  of  motion  but  is  valid  only 
when  the  geometry  is  symmetric  and  the  magnetic  field  is  so 
weak  that  it  does  not  affect  electron  orbits.  The  analytical 
technique  has  been  extended  to  more  general  geometries,  but 


still  without  magnetic  field  effects  Calculation  of  densities 
from  the  equations  of  motion  (e  g.  following  particle  trajec¬ 
tories)  is  valid  for  all  geometries  and  magnetic  fields  but  is 
frequently  hampered  by  numerical  approximations  and  statis¬ 
tical  noise.  Both  techniques  are  used  in  the  calculations  below 
Analytical  expressions  are  used  in  the  bulk  of  the  calculations, 
and  to  check  those  results,  a  few  calculations  were  done  by 
following  electron  trajectories. 

In  general,  the  sphere  collects  only  a  fraction  of  the  electron 
current  which  enters  the  sheath  Even  with  no  magnetic  field, 
electron  angular  momentum  can  prevent  sheath  electrons 
from  reaching  the  sphere.  However,  for  the  46-kV  potential 
and  the  plasma  parameters  being  discussed,  this  is  not  a 
factor;  all  the  390-mA  current  entering  the  7.1-m  space- 
charge-limited  sheath  would  be  collected  in  the  absence  of  a 
magnetic  field.  A  magnetic  field  can  reduce  the  current  col¬ 
lected  Self-consistent,  steady  state  calculations  (Vf  J  Mandell 
et  al„  Current  collection  by  a  high-voltage  sphere  from  a  cold 
magnetoplasma,  submitted  to  Journal  of  Geophysical  Re¬ 
search,  1988)  show  that  high-voltage  spherical  probes  in  iono¬ 
spheric  plasmas  collect  currents  close  to  the  limit  derived  by 
Parker  and  Murphy  [1967],  The  Earth's  magnetic  field  of  0  4 
G  would  reduce  the  collected  current  to  only  14  mA  from  the 
390  mA  corresponding  to  the  space-charge-limited  value.  Ac¬ 
cording  to  recent  estimates  |P.  Palmadesso.  private  communi¬ 
cation.  1988),  turbulence  may  cause  cross-field  transport  and 
increase  the  collected  current,  but  it  would  still  be  only  a  small 
fraction  of  the  sheath  current.  The  calculations  of  the  current 
to  the  spheres  presented  here  all  include  the  effects  of  the 
magnetic  fields  on  the  electron  trajectories.  Jhe  magnetic  field 
is  particularly  important  for  small  rocket  body  potentials  and 
large  electron  sheaths. 

Around  a  negative  sphere  an  identical  sheath  would  form, 
but  with  space  charge  of  the  opposite  sign.  Ion  motion  is 
essentially  unaffected  by  the  Earth's  magnetic  field  When  kil¬ 
ovolt  ions  impact  aluminum  oxide  surfaces,  they  generate  a 
large  number  of  secondary  electrons  [Diet:  and  Sheffield. 
1975].  For  the  high  voltages  experienced  during  the  SPEAR  I 
flight,  each  ion  most  likely  produces  more  than  10  electrons. 
Secondary  electrons  produced  by  ion  impact  dominate  the 
current  to  negative  potential  surfaces.  However,  since  the  elec¬ 
trons  accelerated  by  the  high  fields  travel  more  than  170  times 
as  fast  as  O'  ions,  the  ions  still  dominate  the  space  charge 
density  around  a  negatively  charged  object: 


The  secondaries  are  accelerated  along  field  lines,  and  almost 
all  pick  up  too  much  angular  momentum  to  be  collected  by 
the  spheres.  The  calculations  below  include  the  current  from 
ion-produced  secondary  electrons  but  ignore  their  contri¬ 
bution  to  the  space  charge 

Einally.  when  ihc  spheres  are  biased  positive  with  respect  to 
the  rocket  bod;.,  the  rocket  body  floats  lo  a  negative  potential 
wuh  respect  lo  the  plasma  It  floats  to  the  potential  where  the 
ion  and  secondary  electron  current  to  the  body  exactly  bal¬ 
ances  the  electron  current  collected  by  the  spheres  The  iono¬ 
spheric  plasma  currents  respond  to  a  change  of  body  potential 
on  the  time  scale  for  a  particle  to  transit  the  sheath  less  than 
a  microsecond  for  electrons  and  less  than  a  millisecond  for 
ions  In  comparison,  the  voltage  is  applied  lo  the  spheres  for  a 
second.  The  floating  potential  response  is  so  fast  that  at  least 
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POL  AR  2  *1  plasma  sheath  The  larger  ion-collecting  sheath  partially  engulfs  the  electron-collecting  sheath 


in  a  time-averaged  sense,  the  net  current  to  SPEAR  1  must  be 
zero 

While  the  elements  of  the  theory  presented  here  are  not 
ness  some  components,  particularly  the  mechanisms  control¬ 
ling  electron  collection  from  the  ionosphere,  have  been  at  odds 
with  data  from  lower-altitude  rocket  experiments.  The  calcula¬ 
tions  that  follow  show  that  all  three  elements  (space  charge 
sheaths,  collected  plasma  currents,  and  current  balance)  are 
important  in  order  to  understand  how  the  ionosphere  re¬ 
sponded  to  SPEAR  I 

CaI.CL  1  ATIONS 

The  calculations  presented  below  were  performed  to  com¬ 
pare  the  flight  measurements  with  the  classical  theories  of 
electron  collection  from  a  magnetized  plasma.  The  calcula¬ 
tions  were  done  using  the  NASCAP  LEO  and  POLAR  2.0 
computer  codes  which  take  account  of  the  SPEAR  l  three- 
dimensional  geometry  as  well  as  space  charge  and  magnetic 
fields  While  the  experiment  ideally  would  consist  of  spheres 
floating  in  space,  the  actual  experiment.  Figure  I,  consists  of 
two  spheres  supported  by  meter-long  resistive  bushings 

N  \SCAP  LEO  [  Vf atulcll  a  al.  1982]  and  POLAR  [Cooke 
cl  ■</..  1985.  Lilli- y  cl  ill .  1986.  Katz  cl  al..  1984]  both  solve 
Poisson's  equation  in  three  dimensions  to  determine  the  elec¬ 
tric  potential  distribution  in  the  plasma  surrounding  a  space- 
cralt  They  differ  in  that  POLAR  uses  particle  pushing  to 
sob c  for  the  self-consistent  space  charge  density  by  iteration, 
and  N  ASCAP  LEO  uses  space  charge  obtained  from  formulas 
based  on  current  continuity  in  a  homogeneous  background 
plasma  Both  codes  use  particle  pushing  to  determine  the 
plaiina  currents  Substantial  testing  has  shown  that  in  the 
absence  of  a  magnetic  held,  the  formulas  in  NASCAP  LEO 
have  a  wide  range  of  applicability  and  can  be  as  accurate  as 
the  particle  pushing  results  By  using  the  analytical  formulas. 
VAS(  AP  LEO  can  find  the  sheath  in  less  than  one  tenth  the 
time  of  the  self-consistent  POLAR  code.  However,  for  the 
unusual  geometries  and  strong  magnetic  effects  seen  in 
SPE  AR  I  it  was  necessary  to  perform  self-consistent  POLAR 
calculations  for  at  least  a  few  cases  in  order  to  assess  the 
.ii.cur.icv  of  the  more  approximate  NASCAP  LEO  calcula¬ 
tion' 


Both  codes  use  the  sharp  sheath  edge  approximation,  which 
defines  the  sheath  edge  to  be  the  potential  contour  beyond 
which  the  repelled  species  is  entirely  excluded  The  Debye 
length  for  the  E  region  ionospheric  plasma  is  extremely  short, 
of  the  order  of  a  centimeter,  and  the  sheath  dimensions  are 
la'ge.  of  the  order  of  tens  of  meters.  The  zone  size  at  the 
sheath  edge  is  of  the  order  of  half  a  meter.  To  compensate  for 
the  large  zone  size,  the  sheath  boundary  potential  was  taken 
to  be  10  eV  The  grid  for  both  codes  consisted  of  several 
embedded  grids  of  progressively  finer  resolution.  The  finest 
grid  resolution  is  5  cm  and  is  used  to  resolve  the  space  around 
the  spheres  Several  calculations  of  current  collection  by  a 
sphere  were  performed  using  the  three-dimensional  codes  and 
compared  with  results  from  lower  dimensional  analyses.  The 
two-dimensional  simulations  including  a  magnetic  field  per¬ 
formed  for  this  comparison  are  discussed  elsewhere  (M  J 
Mandell  et  al..  Current  collection  by  a  high-voltage  sphere 
from  a  cold  magnetoplasma,  submitted  to  Journal  of  O eophysi- 
i al  Research.  1988)  The  three-dimensional  results  are  within 
20"..  of  the  comparable  two-dimensional  calculations.  For  the 
cases  in  which  magnetic  fields  are  important.  POLAR  is  in 
good  agreement  with  the  two-dimensional  results,  while 
NASCAP  LEO  makes  errors  as  large  as  50"...  POLAR  is  able 
to  correctly  model  the  contraction  of  the  sheath  perpendicular 
to  magnetic  field  due  to  v  x  B  forces  on  electrons;  the  sheath 
NASCAP  LEO  calculates  is  unmodified  by  magnetic  fields 
The  same  ionospheric  plasma  parameters,  n,  =  5  x  10"’ 
in  \  kT.  ^  0.1  eV.  and  |B|  =  0.4  G.  were  used  in  all  the  calcu¬ 
lations  These  parameters  are  consistent  with  the  flight 
measurements  In  both  codes,  plasma  currents  are  calculated 
by  tracking  test  particles  in  from  the  sheath  edge  until  they  hit 
the  ohiect  or  leave  the  sheath  The  particles  are  subject  to 
both  electrical  and  magnetic  forces.  Owing  to  the  magnetic 
held,  a  small  fraction  of  the  electrons  neither  leave  the  sheath 
nor  arc  collected  by  the  object  even  after  extensive  tracking 
Our  two-dimensional  analysis  has  shown  that  self-consistent 
space  charge  contracts  the  sheath  perpendicular  to  the  mag¬ 
netic  field  and  reduces  the  number  of  apparently  trapped  par¬ 
ticles  and  lhal  these  particles  will  eventually  hit  the  object  In 
the  SPEAR  I  calculations  ihe  number  of  trapped  particles 
accounted  for  no  more  than  20",.  of  the  collected  current 
Neither  counting  them  as  collected  nor  discarding  Iheir  contri- 
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Fig  6  Calculated  ion  and  electron  currents  entering  the  sheath 
as  a  function  of  ground  potential  for  ihe  case  where  46  kV  is  applied 
to  one  sphere 


bulion  changes  the  results  significantly;  in  the  results  below, 
they  are  always  included  in  the  collected  current 

Ions  produce  large  numbers  of  secondary  electrons  from 
aluminum  oxide  Normally  incident  10-keV  oxygen  ions  were 
assumed  to  produce  seven  secondary  electrons  when  they 
impact  aluminum  oxide.  The  most  applicable  data  were  for 
Na  '  on  Al  ,Oj  The  value  for  O  '  was  estimated  from  the  Na  * 
results  of  comparable  ion  velocity  Follow  ing  Dietz  and  Shef¬ 
field  the  secondary  yield  •/  was  assumed  to  vary  linearly  with 
velocity: 

7  =  Of  -  f0)  (6) 

where  r0  is  the  threshold  velocity  of  5.5  *  10*  m  s  ~ 1  and  C  is 
chosen  to  make  the  yield  equal  7  for  a  10-keV  ion  The  inci¬ 
dent  ions  were  assumed  to  hit  isotropically,  which  increased 
the  effective  secondary  yield  by  a  factor  of  2.  The  ion- 
generated  secondary  electron  yield  is  the  largest  uncertainty  in 
the  calculations 

The  calculations  simulate  the  SPFAR  I  discharges  when  a 
single  sphere  was  charged  to  the  highest  voltage  and  the  volt¬ 
age  decayed  owing  to  both  the  plasma  currents  and  the  cur¬ 
rent  through  a  fixed  resistor  The  highest  measured  voltage  on 
SPEAR  I  was  45,320  V,  the  calculations  were  done  with  volt¬ 
ages  of  up  to  46,000  V  Four  different  applied  voltages.  46,  24, 
12.  and  I  kV,  were  chosen  to  represent  the  sphere  voltage 
during  the  discharge  For  each  sphere  voltage  a  senes  of  calcu¬ 
lations  was  performed  varying  the  potential  of  the  rocket  body 
with  respect  to  the  ionospheric  plasma  For  each  of  the  rocket 


Fig  T  Calculated  ion  and  electron  currents  collecied  as  a  function 
of  ground  potential  for  ihe  case  where  46  k  V  is  applied  lo  one  sphere 
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Fig  8a  Potential  contours  calculated  by  NASCAP  LEO  for  Ihe 
case  with  one  sphere  biased  to  46  kV  and  the  spacecraft  ground  al 
—  6  k  V. 


body  potentials,  NASCAP,  LEO  calculated  the  ion  and  elec¬ 
tron  plasma  sheaths  along  with  the  current  to  the  rocket.  The 
equilibrium  rocket  body  floating  potential  was  determined  by 
balancing  the  electron  current  to  the  sphere  with  the  second¬ 
ary  electron-enhanced  ion  current  to  the  rocket  body.  This  set 
of  calculations  thus  yields  both  the  potential  on  the  rocket 
body  and  the  current  to  the  sphere.  Calculations  using  the 
POLAR  code  were  performed  for  46  kV  on  the  sphere  and 
both  -6  kV  and  —8  kV  on  the  body.  These  calculations 


I  ig  M h  Puih  of  .in  electron  in  Ihe  poientials  shown  in  Figure  Ku 
Noie  ih.u  ihe  path  is  dramatically  influenced  by  ihe  presence  of  ihe 
ion-collecting  sheath 
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Currents  collected  js  u  function  of  the  ground  potential  including  the  effects  of  secondaries  and  the  presheath 
enhancement  for  the  case  si  here  46  k  V  is  applied  to  one  sphere 


showed  that  for  these  interacting  bipolar  sheaths.  NASCAP 
LEO.  which  uses  space  charge  formulas,  gives  almost  the  same 
results  as  POLAR,  which  uses  space  charge  found  from  self- 
consistent  particle  tracking.  For  a  rocket  body  potential  of 
-  8  k  V,  the  collected  electron  current  by  the  sphere  is  47.6  mA 
(NASCAP  LEO),  compared  with  53.4  mA  (POLAR).  The  col¬ 
lected  ion  current  by  the  rocket  body  is  3  3  mA  (NASCAP 
LEOi  compared  with  3  I  mA  (POLAR).  When  ion-generated 
secondary  electrons  are  accounted  for,  the  floating  rocket 
body  potential  for  current  balance  from  both  NASCAP  LEO 
and  POLAR  is  -8.3  kV 

RfcSt  Lrs 

The  ion-collecting  sheath  surrounding  the  rocket  body 
dominates  the  smaller  electron-collecting  sheath  around  the 
sphere  As  is  shown  in  Figure  5.  even  though  the  body  poten¬ 
tial  is  only  a  small  fraction  of  the  sphere  potential,  the  much 
larger  sheath  forms  because  the  surface  charge  on  the  5-m' 
rocket  body  is  much  greater  than  the  charge  on  the  0.04-m; 
sphere  This  is  only  due  to  their  relative  sizes  The  effect  of  the 
large  ton-collecting  sheath  is  to  pinch  off  the  electron¬ 
collecting  sheath  and  to  reduce  the  area  through  which  elec¬ 
trons  can  enter  from  the  ionosphere  As  the  rocket  body  be¬ 
comes  more  negative,  the  ion-collecting  sheath  grows  until  it 
completely  isolates  the  electron-collecting  theath  from  the 
ionosphere  Figure  6  shows  the  ionospheric  electron  thermal 
current  entering  the  sheath  as  a  function  of  toe  rocket  body 
potential  when  46  kV  is  applied  to  sphete  For  rocket 
bodv  potentials  more  negative  than  —15  kV.  the  lon- 
collecting  sheath  completely  surrounds  the  positive  sheath, 
preventing  any  electrons  form  entering  the  sheath  This  is  a 
lower  bound  on  the  body  potential  and  doesn't  depend  on 
magnetic  field  effects  on  electrons  nor  on  the  yield  of  lon- 
generjted  secondary  electrons. 


While  the  electron  current  entering  the  sheath  decreases 
monotonically  as  the  body  floats  negative,  the  electron  current 
collected  demonstrates  a  more  complicated  dependence  on  the 
body  potential  (Figure  7|.  For  small  body  potentials,  the  elec¬ 
tron  current  collected  by  the  sphere  actually  increases  as  the 
current  into  the  sheath  decreases.  The  magnetic  field  reduces 
the  collected  electron  current  to  a  small  fraction  of  the  current 
entering  the  sheath  Ffowever.  as  the  body  potential  becomes 
negative,  it  forms  an  ion-collecting  sheath  which  distorts  the 
symmetry  of  the  electron-collecting  sheath.  This  distortion  en¬ 
ables  electrons  entering  the  sheath  to  v  x  B  drift  into  regions 
of  vtrong  electric  field,  from  which  they  can  be  collected 
(Figure  X).  As  the  body  potential  varies  from  0  to  -4  kV.  the 
fraction  of  electrons  collected  increases  faster  than  the  electron 
sheath  area  decreases  For  large  negative  body  potentials, 
most  of  the  current  entering  the  sheath  is  collected,  and  the 
net  electron  current  collected  decreases  because  the  sheath 
area  decreases  The  potential  where  the  collected  electron  cur¬ 
rent  starts  to  decrease  as  the  sheath  decreases  provides  an¬ 
other  bound  on  the  body  floating  potential  This  bound  does 
not  depend  on  (he  yield  of  ion-generated  secondary  electrons 
but  depends  critically  on  the  effects  of  the  magnetic  field  on 
electron  motion 


TABLE  I  Body  Potential  and  Plasma  Current  Versus  Sphere 
Bias 


Sphere  Bias. 
kV 

Body  Potential. 
kV 

Current. 

mA 

+  46 

-X.3 

45  : 

i-  24 

5.3 

36  9 

*  i: 

1 V4 

-  1 

0.6 

1  4 

Values  are  as  calculated  by  NASCAPLHO 
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The  rocket  body  floats  to  a  potential  between  these  two 
hounds,  depending  on  the  secondary  electron  yields  Since  the 
model  for  the  secondary  electron  yields  that  was  presented  in 
the  preceding  section  contains  a  large  uncertainty,  it  is  impor¬ 
tant  to  consider  the  sensitivity  of  the  results  to  the  yields  The 
three-dimensional  calculations  can  be  used  to  determine  body 
potentials  and  currents  as  a  function  of  the  secondary  electron 
yields  For  the  system  with  4b  kV  applied  that  we  have  been 
considering,  a  yield  of  40  is  necessary  to  produce  an  equilibri¬ 
um  body  potential  of  6  kV  with  a  current  of  7K  mA.  while  a 
yield  of  4  corresponds  to  a  potential  of  II  k  \  and  a  current 
of  20  mA.  Thus  while  the  results  are  dependent  on  the  as¬ 
sumed  secondary  electron  yields  extremely  large  changes  in 
these  yields  are  necessary  to  produce  qualitative  changes  in 
the  potentials  and  currents  I  he  three-dimensional  sheath  cal¬ 
culations  that  have  been  performed  can  be  used  with  more 
accurate  yield  data  if  they  become  available  to  remove  this 
uncertainty.  Figure  9  shows  the  current  to  the  bodv  including 
secondary  electrons  and  the  electron  current  to  the  sphere 
Also  plotted  is  the  net  current  to  SPF  \R  I  I  he  zero  in  the  net 
current  determines  the  body  Moating  potential  When  4b  k V  is 
applied  to  the  sphere,  the  calculated  Moating  potential  is  8  4 
kV  with  a  current  of  4'  mA 

I  able  I  shows  the  calculated  floating  potentials  and  sphere 
currents  lor  lour  sphere  bias  voltages  In  Figure  10  the  sphere 
currents  are  compared  with  a  tit  i»v  the  Might  data  \Im.  shown 
is  a  line  through  the  ougin  lit  to  the  four  calculated  points  A 
straight  line  also  tils  the  Might  data  I  he  straight  line  through 
the  calculated  current-voltage  points  provides  an  elleclive  re 
sist;!«*v  '  which  can  he  compared  with  the  same  quantity  from 
the  Might  data  The  apparent  linearity  <4  the  current-voltage 
relationship  is  peculiar  to  the  SP1  \R  I  conhguiation  and  the 
voltage  range  of  the  calculations  The  agreement  hetween  the 
calculations  and  experiment  is  well  within  the  estimated  un¬ 
certainty  of  the  calculations,  which  is  about  40  Figure  4 
shows  the  calculated  Moating  potentials  compare!  with 
bounds  determined  from  ion  spectra)  measurements  made 
during  Might  As  discussed  earlier  the  spectra  were  much 
broader  than  anticipated,  possibly  as  a  result  of  backscatter 
from  (he  rocket  body  or  from  sheath  oscillations  The  calcula¬ 
tions  all  fall  within  the  experimental  hounds 


Disci  ssio\ 

SPLAR  I  measured  the  response  of  the  ionospheric  plasma 
to  l.trge  obiccts  with  exposed  high  voltages  The  rocket  Mew 
high  enough  and  released  little  enough  neutral  gas  that  ioniza¬ 
tion  within  the  sheath  did  not  cause  the  high-current,  low- 
voltugc  glow  discharges  observed  on  low -altitude  electron 
beam  rockets  and  also  observed  during  large  chamber  testing 
ol  SIM  \R  I  The  calculated  Moating  potentials  and  currents 
collected  are  in  agreement  with  potentials  and  currents  mea- 
siired  during  the  Might  That  the  most  important  physical  pro¬ 
cesses  are  included  in  the  calculations  is  supported  by  the 
close  agreement  hetween  the  calculations  and  the  experiment 
in  such  a  complicated  system  High -impedance  space  charge 
sheaths  were  maintained  until  the  rocket  passed  below  !  40  km 
ui  altitude,  where  ionization  effects  become  dominant 
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The  POLAR  Code  Wake  Model: 
Comparison  With  in  Situ  Observations 


(  i  Ml  KI'H  1 

. . .  /  .if’ . .  I 

I  Km/ 


V(  l«.  ■/  /  1  (  I. 


Mc.isar emeu -  ol  :hc  .on  .aid  eles Hi m  dens iiic's  i.m\ lateit  .\ i! h  'ho  -v  Ac  - '(  he  -him k  oi  tali .  .uu 
niaje  Kv  'he  Plasma  Diagnostics  Pa.k.ige  '  PDP'  limine  ihc  !‘iX'  Sp.uelaP  2  ni"i'in  i  ;.<s-  -t-eion-.  .  I 
■  he  *ake  ..I  ilio.nices  ol  MS  2^0  m  Ji>\\n>irc.im  .irul  nic.i^iiuMnom^  .long 'he  a. me  i\i.  Irani  '  * . »  'mi 
ill  were  .’Plained  I  ht*  Pill  VR  *akc'  inode!  devc-tnped  loi  Ihc  \u  foi.e  (icopni-ii-  1  ipoi  .loiv  :*< 
.tudv  .h.ngine  ■'!  -paiL-erall  in  low  ilnit.de  high  I  in.ii  i.  .n  .'ifih  was  u-cil  V  Pcitoii'i  i  'hi  ..-c 
dimensional  'in.u.aiion  .  M  Ihc  plasma  wake  evaluated  .if  points  .lone  -cl. .Inc  uceori  .'1  he  PI  >P 
The  POI.  \R  . i <<Je  me'  several  simplilving  assunipnons  »  picjic!  wake  den'."es  PiC'-c  ni!nde 
negle*iine  ih.  pv.cneiu  held  .mil  assuming  'h.ti  'he  plasma  .s  e^u. iM-neuif.il  I  he  .ode  noilels  plasma 
deti'Uv  .the. id  el  Ihe  expansion  tiont  us.ne  i  ncati.il  appiov.m.uion  .old  mode's  '.he  pi. mn. .  .lensiii 
behind  ’he  esp-awon  lieni  K.  imne  'he  .ell  -mnlai  solniion  >•!  'he  expansion  .*!  i  pl.i.m.i  nit .  i 
vacuum  for  .uses  -s he i e  /  /  ihe  Hunt  .s  mu  .h.u r  .uni 'hci  nial  motion  . an  u.otmliot  mosi  id 

Ihe  expansion  I  hi.  ipproa.n  is  .ompulahon.illv  .civ  ctli.icni  I  he  iC'iills  presented  he'e  .tie  ihe  hi  si 
known  ...nir,, risen  k'ef.eee*i  -u.h  i  m-ilc!  ifid  aifrl.1  ft  ..ru  d.d.i  .mr. tilled  lei  ■  'hie.’'  el  s..,ic  -i/e 
in*  '  I  t.eli'f  gii.iiii..'.u'  I'tdji.i'tiii.i'.'c  iccietk  f  .icH-uriil  ri:  ilmau.c  c'..'i*  h  .n  ■"  .'i 
:  nilu.d  i  ne  ih.-.i  .d  :c.m  'e  * '  i  - 1  elder  me  neilc!  ..rp  .mii!..'.i.  *n.  .,i  e  .el  ilie.l  \:i  ,ni:  emi:  .1.  -r  a  u  - 
between  'he  nedci  mil  •!  .1..  sucgesis  ih.d  'he  oifiie:  •  ne  a  -'.ike  max  ne  hik'd  hied. 'Him. .mix  hx  , 
pick-up  nn  piipie  ihun  .reeled  hem  neuitai  .om.immanis  .nil  ’hei  *he'C  x-uiid  '  ne  ’e  -c  n.'i.ded  i! 
.iei.iir.ile  w.,Ce  models  el  .;fee  gav-cnmtine  .'I'iClIs  ere  leqrmed 


Is  I  K<  >IH  t  r  ids 

In  this  paper  ne  discuss  measurements  made  bv  ihe 
Pl.ismj  Diagnxisiics  Package  i  PDF'  during  Spacelub  2.  which 
.ire  presented  bv  \f,irph\  cl  ,ii  1 1 s»X*-* | .  ,rnd  Lompure  ihose 
rc'dlts  with  predictions  triirn  (he  Air  Force  Geophwics 
l.jhor.ilorc  i.AFGI.i  POLAR  w.ike  code,  which  uses  a 
.ornplec  geometric  model  ot  ihe  oibiter  and  ihe  selt-similar 
solution  ot  (he  expansion  ol  a  plasma  into  a  vacuum  as  ils 
model  basis  Previous  reports  | K,il:  cl  til.  IW.'I  have 
.ompared  the  predictions  ot  Pt  tt.AR  lo  observations  ot  T 
I  plasm.is  in  ihe  laboraiorv 

Caution  shv'uld  be  exercised  in  extending  conclusions 
iboui  ihe  iccuracv  ill  ihe  POI  VR  model  lo  eonclusions 
regarding  veniicurn*"  ot  the  underlying  phvsical  processes  ii 
contains  Several  other  invesiiganons  nave  studied  (he  ap- 
plicubilitv  ot  the  self-similar  malhem.iiic'  |.Viimi>  ci  ,il  .  I9XV 
Nil  w  haiijliiin  cl  nl  .  IdSh.  I i,n c<  it /i  d  ill. .  Idh9:  Du  hi’lil  a 
ji  .  Ids'  Kn:ima  cl  ill ..  ldHS|  (o  wakes  Ii  is  our  purpose 
only  u>  determine  it  Pt  >L.\R  provides  a  reasonable  mi'del  lor 
the  wakes  ot  "large  ‘  obieets  in  the  ionosphere,  as  it  has  tor 
I  I  plasmas  in  the  laboraiorv 

We  describe  brietlv  the  POLAR  model  and  review  the 
phvsics  it  contains,  .ompare  the  data  wi(h  the  model,  and 
■hen  discuss  the  range  of  validity  of  the  code. 

Tm  P< iL.‘.R  Loot 

To  develop  a  code  that  can  adequately  describe  the  plasma 
wake  behind  a  large  object,  particularly  one  of  complex  gcom- 

t.i’pvright  IdHdbv  ihe  Vmenc.in  (jeophvsieal  f  mon 

Piper  number  S4J  -MMIXS: 
max  .ii”*  xu  s,JJ  A-uo.''X_"yi*  iki 


ctrv  i arcful  consideration  ol  assumptions  and  approximations 
arc  required,  as  arc  simplifications  allowing  tor  computational 
efficiency  Ihe  POLAR  code  has  evolved  with  such  lon'ider- 
alums  in  mind  \  detailed  description  ot  Ihe  POLAR  wake 
model  is  given  by  Kit/:  of  <il  j  |UR'|.  and  it  is  ihe  purpose  here 
onlv  lo  review  the  basic  phvsie.  and  processes  in  POI  \K  so 
the  reader  m.iv  have  some  insighl  into  the  vaiiditv  ol  the  eode 
I  he  model  of  the  wake  structure  used  bv  POL  VK  depends 
on  ihe  position  relative  lo  ihe  so-called  ion  tronl  I  his  ion 
from  marks  the  bound. irv  where  electron  Jensitv  begins  to 
change  on  a  scale  commensurate  with  the  Dcbve  length  and 
the  ion  densitv  takes  a  sudden  and  dramatic  drop  Several 
authors  have  discussed  Ihe  relationship  bgt w ecu  ihe  wake  Mil 
process  and  ihe  thcorctiv.il  piohlem  ol  ihe  expansion  ol  a 
plasma  into  a  vacuum  In  paitiiular.  problems  applicable  to 
ionospheric  conditions  have  been  Heated  bv  <i,n on  i  r  ,u 
JIMbfij.  lltncMiii  iiiitl  I'ihh  '.  \k;i  |  I4's|.  md  So’e'i  ami 
Si  hunk  |14«2|.  to  name  a  lew 

The  solution  te>  the  Vlasvciv -Poisson  equation  svstem  is  in 
general  quile  ditficull  lo  obtain,  bill  lor  the  expansion  ot  a 
plasm. i  into  the  void  it  l  an  be  sob  ed  e  xplicillv  1  (  nu .  ■ n  h  ,  i 
til..  IUbd|  Ahead  of  Ihe  ion  lioni  the  plasma  is  Heated  as 
rarefied:  its  motion  is  c<>nlrolled  bv  the  thermal  spread  in  ion 
velocities  Behind  the  front  the  motion  is  eonnoiled  bv  ihe 
electron  temperature  and  ion  m.:"  I  icurc  I  illustrates  lhc'0 
regimes  and  defines  ihe  eooidin.iie  svstem  used 

The  governing  equations  in  the  region  behind  Ihe  fiont 
considering  th.it  clcelrons  are  more  mobile  than  ions  and  that 
thei  maintain  equilibrium  wuh  a  local  potential,  are 

The  Bolt/man  relation 

exp  <  fAV  /  i  ill 
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* POLAR"  WAKE  MODEL 
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Fig  I  Inc  K)L  AK  wake  ctnlc  disltnguishes  three  regions  of  interest  Ihe  ambient  plasma  the  region  of 
self-similar  model  and  the  neutral  approximation  spaces  are  hounded  hv  the  Mach  cone  Z  -  -  S.,i  and  ion  front 
respectively  The  coordinate  s\stem  used  is  consistent  wiih  equations  ill  through  <l(h 


C  onlinuttv 


[qua  l  ion  ot  motion 


Poisson'  equation 


rUl,  (Hll.t  I 


;  n;  -  ftut> 

Tr  '  “  ‘  vTT 


4  Tl  I  II  -It  1 


Km:  ft  til  |  l-4Ns]  showed  lhal  this  furmulj  agree'  well 
with  laboratory  data  Irom  W'lvlu  rr  ill  [l-4X'j  and  incorpo¬ 
rated  it  in  POLAR  *head  of  l hi"  front  /,  the  pl.i'tn.i  i' 
assumed  to  expand  owing  to  thermal  motion,  the  so-called 
'  neutral  approximation  Behind  /,  the  plasma  evoKc' 
into  a  state  which  is  self-similar  |  Chun  it  ill  .  I -JS4 ,  I  ne 
self-similar  solution  of  i  |c--4i  toi  :  ■  V, ,/  is 


where  S . ,  -  tkf  Mr  '  is  the  ion  acou'-.i.  speed 
-■*'  Ihe  lime  xaiiahle  o  dehned  .o 


n,,  jmhiem  densitx . 
n  ion  densitx  . 
n,  electron  densitx 
/  electron  temperature 
i  electron  charge 
it<  loc.d  potential . 
k  Bolt/man's  constant 

and  where  .  is  a  xariable  representing  distance  parallel  lo  ihe 
front  velocity  or.  in  this  case  perpendicular  m  ihe  orbital 
velociix 

(r,m  ft  ill  |  |M?5|  have  numerically  solved  1 1 1  through  i4i 
to  predict  the  position  of  the  ton  Iront  Km:  <i  •  il  IIW'I 
developed  an  analytical  ht  to  the  (  row  results 


-  V 


f  ■*"  '  L  ) 


I  n  i  I  -  t  r  utf 1  (Of 


0  4  Nw  I 

I  —  —  1 1  itj/  Ini*  mi//1 
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xx here  '  is  the  distance  behind  (he  oh icsi  iperpcndicul.ir  U»  ■ 
sind  \  ,  is  ihe  orbital  \cioot\  He  define  ihe  self  simii.o 
% ariahle  _  j" 


I  hit"  ihe  self  similar  solution  cssentialk  stales  that  he 
(ween  ihe  region  hounded  in  positive  :  h\  (he  from  / ,  and  in 
negative  :  hv  ihe  line  S ihe  densttv  rises  exponent- 

tiall\  lo  he  equal  lo  (he  amhient  value  along  ;  -  -  S./  I  hi" 
is  an  iniuilivek  re.isonahle  result 

In  summary,  ihe  wake  routines  in  Pl)l  emplov  two 
limnini:  cases  <!i  Ahead  of  ihe  ion  Iron!  ihe  eleefriv  held  i» 
negligible  and  (he  motion  of  ions  jv  ideniie.il  fo  neutrals  O 
behind  ihe  ion  Iron!  whose  position  is  determined  hv 
the  qiixisi-neuK a!  sell-sinulat  solution  of  <bi  is  implemented 
POl  \K  has  routines  whuh  model  ue v tiraiclv  the  geome 
I  r  \  yd  (he  uhiut.  .mi1  *he  tuva!  u»n  1 1  .ihl  u 11  ie '  "e 
saL  111. tied  Ill'll'. 


are  the  ion  plasma  frequent. \  and  |Xrs\e  length  respev 
(ivek  and  h  is  a  lice  paiamelef  dcici  milled  lo  he  I  b 


where  ‘ is  ifu.  unpvifuihed  JM '  i^uiion  t  n  :  i«  »n  to:  ,t 
dnluni:  Maxwellian  and  c  \  i  2 1  has  value  »f  *'  lu 

starling  ftoni  \  aiul  eoine  in  ihe  dneslion  l.1  would  si:  ike  ihe 
v  ehu  le  anel  I  if  ll  v\  t  *uld  u*  >1 
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X| ,  arm  K  "Z  i*‘ 


iM  t  . 


VV 1 1 1 1  1  ' ' ‘ - 1 


f  t  '  l  xirtk!  'he  relative  ’raicwtorv  al  ’he  CUP  and  arbiter 
i>(  t  \K  *aiviiialcs  the  densi'v  as  a  turu'ion  id  time  toi  'he  ennie 
’•ee  ’'  em  '.me  penou  \re.o  •  *1  m'cresi  jre  'aheleJ  in  the  rigure  !  he 
1-  ."  eJ -'arm. ili/.mon  values  for 'he  pl.i.m,i  arc  n  I  •  III  'em 
1  ;;IKI  K 


’  he  'ocal  densitv  is  given  h\ 


I  lu.  milial  densitv  calculated  in  three  dimensions  lor 
neutral  particles  in  compared  with  density  calculated  .innuiu- 
me  ’he  complex  geometric  object  in  repl.e  cd  by  a  flat  plate  at 
a  position  where  the  dominant  source  appearN  at  the  object 
edge  This  ratio  provides  a  "geometric  correction  factor, 
■which  in  applied  to  the  quasi-neutral  one-dimensional  solu- 
tion  discussed  earlier  for  positions  behind  /,  In  this  wav. 
I’t  >1  \K  can  calculate  e)uite  rapidly  an  approximate  value  tor 
tne  ion  and  electron  densities  in  the  wakes  of  complex 
objects 

'sole  that  the  assumptions  behind  the  front  are  1 1 1  that  the 
electron  temperature  and  ion  mass  govern  ihe  equation  of 
motion  1 2 »  that  the  plasma  is  quasi-neutral.  1 3 >  that  the 
magnetic  held  does  not  ailed  the  ion  or  electron  motion,  t4l 
that  equation  i'i  serves  as  a  cov'd  approximation  for  deter¬ 
mining  the  boundary  of  the  ion  front,  and  1 5 »  that  the 
geometric  correction  factor  calculated  in  detail  with  the 
’lircc  dimensional  neutral  model  can  be  approximately  ap¬ 
plied  to  correct  the  plasma  densities  as  well  Therefore  the 
neoruhm  c.m  address  complex  geometries  but  takes  advan- 
'.ice  at  the  smooth  wake  structure  characteristic  of  iono¬ 
spheric  plasmas  where  T  T  I  \ddilionally .  the  mvidel 
’Cttores  fields  existing  in  a  sheath  nejr  the  body  surface, 
w  nic h  'hould  not  he  of  concern  in  our  case  w here  the  orbiter 
is  near  plasma  potential  This  implies  that  ion  acceleration 
-  iic  il, iieil  hi  P<)|  AR  is  dominated  bv  electric  fields  due  to 
'P  ice  ch  iree  separation  in  the  wake 


Gjxip  xrisos  or  R i  s l  l  is 

I  cure  2  is  a  plot  of  the  normalized  densitv  as  a  function  of 
'lire  predicted  bv  POI.AR  for  the  PDP  orbiter  relative 
traicclorv  The  density  ratio  V„  \  ilnhlvl)I  varies  from 
aprroximaielv  DR  for  wake  transit  iWTl  I.  which  is  at  a 
distance  of  24s  m  to  ii  s  jlir  t p,e  v. lose-st .  W  f  }.  at  4S  m 
1  tne  musi  be  careful  when  comparing  the  model  to  the 
data  since  the  model  assumes  a  tixed  background  densitv  ol 


III  whereas  '.tic  ictnal  loiio'pnciic  •Ifti'if.  c  m  m.l 

ill  i  c  s  ,||V  v  oil  s  Hi  c  1  a  I '  I  V  I  !  1  c  I'ac  k  C '  OH  lul  Jc  II "  1 1  i  ':,"01  |,  'i 

:hc  lllovlel  IS  ".('leal  ol  llial  oi’scl'-cil  III  till  :  o.  nil  . . 

spticic  1  tic  I'salUCil  model  tcilpci. lime  IS  .  "  lalll 

2'im  K  Ih-:  otscivc'd  'c  mpci  almc  vane's  pno  ->i  uinu, 

2 ' ‘  ,  1 1  om  t  tio  assumed  v  .due  dut  me  t idles  at  uitc  t  c  -I  !  :ie 
wake  transits  w e i c  all  planned  io  occiu  in  the  ret.nivelx 
stable  davside  umosphere.  and  over  the  time  'P.m  ol  i  ci'tn 
wake  transit  the  background  densitv  is  believed  to  be  stable 
to  -  Ilf  ;  lUurpIn  (7  ill  .  IRHR|  In  the  following  discussion 
we  will  alwavs  cv'iupare  the  observed  densitv  ratio  durinc  i 
w.Ao  tiansit  to  the  percentage  chance  piedictcd  bv  the 
mi  >del 

I  iciiic  'ii  plots  ihe  I'h.crved  electron  anil  ion  densitv 
depression  during  the  hack  away  maneuver,  with  a  constant 
It)'' cm  as  i he  leleience  densitv  |  \hnpli\  imi  .  1'i.Vt!  I  he 
densitv  chosen  lor  normalization  is  atbiltaty  and  does  not 
atiect  the  shape  of  the  profile  as  long  as  It  IS  constant  ll  Is  Ihe 
shape  and  relative  depth  of  the  profile  we  wish  in  compare 
with  P<  >LAR  I  he  dots  in  Figure  }<i  are  the  electron  densitv 
ismall  periodic  depressions  result  when  the  boom- mounted 
probe  passes  through  the  wake  of  the  .pinning  PDPi  Ihe 
plus  svmbols  are  the  ton  densitv  obtained  once  per  spin 
cycle.  During  the  hack -aw  ay .  the  PDP  and  nrbucr  pass  if,, m 
a  latitude  ol  approximate!'.  -  2“  to  ■  '-(>  meailx  altemoon 
local  tunc  I  icurc  'a  also  shows  ihe  piedicled  noim.ilizcd 
Jcnsiix  calculated  bx  the  Rl)|  \R  model  described  caiher 
\s  can  be  seen  good  agreement  exi't'  at  distances  liunt  -ll 
to  -  *'  m  It  should  also  be  noted  that  ovet  the  range  ot  the 
data  shown,  the  RPA  and  I.P  data  agree  on  densitv  .vithin 
-  MV  , , 

As  discussed  bv  l Unpin  or  til.  [|RSR|.  we  cannot  ea'ilv 
normal.ze  to  either  the  prev  ious  or  follow  me  orbit  so  we  U'C 

a  constant  value  of  HD  cm  'tor  V  . . .  This  .alue  ippeai' 

to  be  the  peak  densitv  observed  alter  the  PDP  has  Ivlt  ihe 
wake  approximated  I  hour  and  40  mm  later  Over  the 
relatively  long  period  of  this  maneuver  the  inno'phe'e  ".as 
siirelv  changed  bv  more  than  l(f  <  and  this  mav  explain  'he 
disparity  alter  /  20  min.  since  the  PDP  to  otoiict  distance 

is  changing  relatively  little  figure  'b  is  a  pioi  .|  : he 
predicted  densitv  calculated  bv  the  IRI  lonospheic  node! 
and  the  measured  densitv  one  orbit  later  'same  oc.u  '.me 
and  latitude'  Note  that  the  gradual  satiations  observed 
owing  Io  changing  local  time  and  latitude  are  consistctil  w  1 1 U 
Ihe  general  trend  predicted  bv  the  IRI  model 

Let  us  turn  now  to  the  wake  transit  observations  I  able  I 
lists  ihe  wake  transits  md  compares  ihe  prediv'cd  and 
observed  depletions  as  weil  as  conditions  ,u  the  time  ot  the 
center  of  the  wake  Note  that  in  this  enure  r  mce  ol  'iu2'0 
m.  the  calculations  and  observations  agree  to  within  ;(f  ■ 

Since  the  orbiter  has  a  complex  ecometrv  the  detail'  ol 
the  wake  structure  at  a  distance  of  4' m  iW  I  ' 1  slow  mticam 
may  be  affected  Figure  4  plots  a  detail  of  that  vv.,kc  'r.insit 
and  the  POL  \R  model  as  a  function  ot  tunc  I  he  data  ate 
normalized  to  a  constant  value  since  it  is  clear  some 
variation  in  background  densitv  occurs  bul  it  is  not  dear 
exactly  what  lhal  variation  is.  I  o  normalize  to  am  unknown 
value  other  than  a  constant  would  introduce  artificial  varia¬ 
tion  and  skew  the  results  Ump/n  i  t  at.  |  |R«4|  discuss  this 
issue  in  detail  and  in  lhal  case.  XV  1  ’>  is  normalized  bv  the 
data  from  the  prior  orbit  However  \1m pin  i  i  til  |l’<XR|are 
attempt  me  no  comparison  to  a  model  Dunne  this  wake 
transit  I  he  bac  kgiound  is  be  I  lev  ed  Io  cm  bv  as  much  as  Ilf. 
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TIME  /MINUTES 
GMT  HOUR  0 


Fig  }  Ion  and  electron  data  during  the  back-away  are  ploiled  in  Figure  V.  normalized  to  «i.  *  I  •  llT  cm  Sole 
the  relatively  good  agreement  between  model  and  data  between  '0  and  m  Heyond  r  -  ft)  nnn  i'*  mi  background 
density  varies  considerably,  as  is  illustrated  in  Figure  .*/>.  which  include'  data  tor  one  orbit  later  than  Figure  so  at 
approximately  the  same  station-keeping  position,  superimposed  on  a  predicted  densitv  tor  this  orbit  from  the  IRI 
ionosphere  model 


| Murpliy  cl  ul  .  I989|  so  the  model  cannot  be  tested  to  an 
accuracy  greater  than  that. 

Discussion 

In  examining  the  back-away  densitv  profile,  we  find  three 
relevant  observations  from  Figure  V 

I  Close  io  the  orbiter  (<  TO  ml  the  model  underestimates 
the  observed  density  by  I  to  2  orders  of  magnitude. 

2.  In  the  range  30—75  m  the  model  predicts  quite  accu¬ 
rately  the  gradual  increase  in  densitv  until  the  time  /  20 

nun  (Figure  3).  Alter  /  -  20  nun  the  observed  densitv  seems 
to  have  a  variation,  which  is  not  believed  to  be  wake-related 
I  hese  densitv  changes  result  I  torn  lono'phciic  variability  as 
the  spacecraft  approaches  the  daw  n-dusk  meridian  plane  and 
are  predicted  bv  empirical  models  such  as  I K I 

Considering  the  first  observation,  recall  that  the  assump¬ 
tions  incorporated  within  the  POI.AR  wake  model  req "ire  a 
quasi-neutral  plasm. i.  assume  a  sell-sinnlat  solution  and 


neglect  magnetic  fields  Since  the  electron  and  ion  densities 
observed  even  at  the  beginning  of  the  release  and  back-away 
period  agree  within  I0'i.  it  would  seem  that  quasi-neutrality 
would  be  valid  It  has  been  shown  by  Chan  cl  ill  |l984|that 
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•  ^angmlhr  probe  data 

♦  POLAR  model 


14  i-  14  i'  1 4  19  50  5'  ^  J  53  5i 

GMT  4h  '  ME.  '1  vt''.  .  -ii 

Fig.  -1  Detail  of  data  from  WT  '  is  compared  lo  POLAR  predictions  Data  are  normalized  to  a  constant  Since  there 
ma>  he  variation  in  background  of  10".  differences  ot  th.it  order  between  data  and  model  are  not  signilicunl 


case,  would  always  act  10  limit  plasmTi  flow,  rather  than 
enhance  it  Therefore  it  too  can  he  eliminated  from  serstnu 
as  an  explanation  for  the  poor  tit  at  less  than  .»<)  m 

The  answer  to  the  disparity  between  model  and  observa¬ 
tion  would  seem  to  lie  in  the  role  played  by  contaminant 
tons.  Vfurphv  et  ul.  |  I989|  discuss  these  ions  and  offers  their 
presence  as  an  explanation  for  the  disparity  between  these 
data  and  that  taken  at  similar  distances  while  on  the  remote 
manipulator  system  IRMS). 

Let  us  turn  now  to  the  30-  to  75-m  distance  region.  The 
most  dominant  characteristic  of  both  the  data  and  model  is 
the  relatively  smooth  increase  in  density  as  the  PDP  moves 
axially  along  the  orbiter  wake 

This  midwake  region  has  been  studied  extensively  in  the 
laboratory  and  the  wake-fill  process  depends  strongly  on  the 
body  size,  body  potential,  and  ratio  of  ambient  ion  to 
electron  temperature  Stone  [I98l|.  Fournier  and  Piutu  he 
[lg~5|.  Hester  and  Sonin  [  I ^''O | .  and  many  others  have 
performed  laboratory  experiments  and  observed  tine  struc¬ 
ture  in  wakes,  including  ion  density  peaks  along  the  wake 
axis  and  wavelike  condensation  disturbances.  It  is  impor¬ 
tant.  however,  to  note  that  for  the  case  of  large  bodies  in  low 
earth  orbit  (LEO)  1  It  the  body  potential  is  not  too  different 
from  the  plasma  potential  la  few  kT  at  most,  since  the  body 
surface  is  an  insulator  and  does  not  expose  v  x  B  potential  to 
the  plasmal.  1 2 )  the  plasma  is  a  cold  Maxwellian  (  lev)  and 
collisionless,  and  <?l  the  ambient  ion  and  electron  tempera¬ 
tures  are  close  to  being  equal 

\n  excellent  review  of  laboratory  work  before  1975  is 
given  by  Fournier  and  Pntaehe  [ I975J.  Another  excellent 
review  of  the  subject  of  expansion  of  plasma  into  the  wake  is 
given  by  Samir  el  at.  ( 1 983 ] .  In  these  cases  the  authors  agree 
with  the  basic  finding  of  Gurevich  and  Pitaevskii  ( 1969)  that 
the  tine  structure  and  ion  peaks  observed  in  certain  labora- 
torv  investigations  vanish  as  T,  approaches  Tr.  We  see  in  this 
case  that  in  spite  of  (or  perhaps  because  of)  the  effect  of 
contaminant  ions,  we  have  a  large-scale  wake  which  is 
basically  devoid  of  any  fine  structure,  at  least  in  the  sense  of 
total  electron  or  ion  density  It  should  be  emphatically 


noted,  however.  1h.1l  this  does  not  imply  knowing  all  there  is 
to  know  about  the  wake  structure  The  overall  plasma 
density  is  onlv  the  zeroth  order  parameter  Ion  and  electron 
composition  \Hrchon  \kv  et  ill..  Ig87|.  vector  measurements 
of  ion  velocities  {Stone  el  at..  1983.  19881:  electron  temper¬ 
ature  | Murphy  et  ill..  1986:  Rain  et  ul..  I984|;  plasma 
turbulence  {Rant  er  at..  I *484)  all  play  roles  in  understanding 
the  total  physics  of  the  wake  structure  for  such  a  complex, 
gas-emitting,  large  object. 

As  discussed  earlier,  the  structural  differences  between 
predictions  and  observations  after  t  ~  20  min  in  Figure  3<t 
are  attributed  to  natural  ionospheric  variation  inot  models  I 
by  POLARl  Figure  3 h  illustrates  the  density  profile  one 
orbit  after  the  back-away  maneuver  and  shows  this  similar 
structure 

Let  us  compare  the  model  predictions  to  observations  for 
WT  I  through  WT  4.  The  agreement  is  qutte  remarkable  and 
affirms  that  the  "well-behaved"  wake  structures  associated 
with  Tt  T ,  =  I  plasmas  can  be  adequately  modeled  by  the 
physics  contained  in  the  POLAR  model.  There  is  onlv  one 
significant  difference  between  model  and  data.  WT  2.  which 
occurs  at  “  125  m.  seems  to  be  considerably  deeper  man  WT 
4.  which  occurs  at  a  little  more  than  100  m  Murphv  et  al. 
1 1989]  discuss  this  extensively,  and  we  do  not  believe  that, 
considering  approximations  made  by  the  model  and  errors 
made  in  normalization,  we  could  expect  any  better  agree¬ 
ment.  If  the  magnetic  field,  contaminant  ions,  and  orbiter 
sheath  do  play  some  role,  it  is  clear  from  both  the  model  and 
the  data  that  it  must  be  a  secondary  one 

Studying  the  detail  of  WT  3  observations  and  POLAR  s 
predicted  profile,  we  also  find  good  agreement.  This  is 
significant  because  at  45  m  downstream  the  details  of  ihe 
orbiter  geometry  and  its  effect  on  the  wake  can  not  yet  be 
"washed  out"  ithe  long  dimension  of  the  orbiter  is  -36  ml. 
The  agreement  between  model  and  data  seems  to  imply  that 
the  geometric  assumptions  are  valid  and  that  it  is  permissible 
to  use  the  geometric  correction  factor  calculated  from  the 
neutral  flow  model,  at  least  to  first  order. 

Note  that  the  center  of  the  predicted  wake  seems  to  be 
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offset  slightly  front  (hut  observed  and  that  the  predicted 
density  gradient  seems  slight! v  greater  than  observed.  Mur¬ 
phy  el  ul.  (1989)  discuss  the  accuracy  with  which  the 
trajectory  reconstruction  takes  place  This  offset  error  is 
consistent  with  that  level  of  precision.  It  is  also  possible  that 
errors  in  normalization  (constant  background  assumed) 
could  produce  this  effect. 

The  difference  in  density  gradients  may  be  due  to  a  slightly 
different  plasma  temperature  than  that  modeled  or  may  be 
consistent  with  the  role  played  by  contaminant  ions  in  the 
neutralization  of  the  space  charge  electric  field. 

Electric  field  data  is  difficult  to  discern  from  the  PDP 
measurements  because  of  interference  from  another  instru¬ 
ment.  but  J.  Steinberg  (private  communication.  March  1988) 
has  examined  data  from  the  time  period  of  WT  3  and  finds  an 
electric  field  which  changes  sign  at  (he  wake  center  This 
field  is  w  uhin  a  factor  of  2  of  that  expected  from  a  self-similar 
expansion.  No  attempt  has  been  made  to  compare  the 
predicted  field  computed  by  POLAR  to  the  actual  data,  since 
error  bars  on  the  data  are  rather  large  The  authors  may 
examine  this  as  well  as  the  RMS  data  \  Trihhle  cl  ul..  1989)  in 
more  detail  in  a  future  paper  The  RMS  roll  data  may  be 
more  useful  for  comparison,  since  those  data  are  taken  at 
<  10  m  downstream,  where  the  density  depletion  and  electric 
fields  are  greater. 

Conclusions 

The  sampling  of  PDP  ion  and  electron  densities  verify  that 
for  ionospheric  plasma  conditions  ( T,  ~  T, ).  the  orbiter  wake 
is  relatively  smooth  in  its  structure  from  =30  m  to  distances 
of  =250  m  downstream.  POLAR  predicts  this  smooth  wake 
structure  and  agrees  with  the  observations  to  an  accuracy  of 
<300  (Table  I)  For  large  and  complex  sy  stems  such  as  the 
orbiler.  outgassed  products  may  play  a  significant  role  in  the 
structure  of  its  wake  at  distances  less  than  ihe  characteristic 
bodv  dimension  Adequate  modeling  in  (his  regime  requires 
input  (hal  details  the  outgassed  species  and  rales,  chemistry 
describing  their  interaction  with  the  ionosphere,  and  inclu¬ 
sion  of  magnetic  field  to  account  for  (he  pick-up  ion  popu¬ 
lation  in  the  wake. 

At  midw'ake  distances  (he  magnetic  field  effects,  if  any. 
would  appear  to  be  secondary  to  the  dominant  role  of  (he 
electric  field  and  thermal  motion  already  modeled  by  PO¬ 
LAR  in  the  wake-fill  process.  However,  lor  slightly  larger 
objects  or  high-inclination  orbits  the  magnetic  field  may  have 
to  be  considered  if  accurate  results  are  required  In  addition, 
the  fundamental  scientific  questions  associated  with  large- 
body  wakes  will  eventually  require  its  inclusion 

For  comparisons  between  in  situ  observations  and  models 
such  as  POLAR  to  be  meaningful  at  levels  better  than  a  lew 
lens  of  percent,  more  than  simple  axial  or  planar  profiles  of 
density  are  needed.  Future  experiments  will  also  require 
good  background  measurements  and  inclusion  of  vector  ion 
velocity,  electric  held,  and  particle  distribution  functions 
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Abstract 

When  subject  to  naturally  occurring  fluxes  of  electrons  with  tens  of  kilovolts  of 
energy,  spacecraft  surfaces  accumulate  negative  charge.  This  charging  leads 
to  potentials  as  high  as  20,000  volts  on  the  spacecraft.  The  potential  achieved 
depends  on  the  electron  spectrum;  the  density  and  energy  spectrum  of  the 
ambient  ions;  the  solar  flux;  and  the  materials,  shape,  and  size  of  the  spacecraft. 
While  the  first  observations  of  charging  were  made  on  satellites  in 
geosynchronous  orbit,  charging  has  been  observed  more  recently  on  satellites 
at  lower  altitudes  in  polar  orbit.  The  physical  mechanisms  which  cause  the 
charging  are  similar  in  both  cases.  Spacecraft  surface  potentials  rise  negatively 
until  the  electron  accumulation  rate  is  diminished  and  the  ion  collection  rate 
enhanced  so  that  the  net  current  to  the  surfaces  is  zero.  In  low  polar  orbit,  the 
enhanced  ion  collection  dominates;  in  geosynchronous  orbit,  both  processes 
are  important. 
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Introduction 

Spacecraft  in  geosynchronous  orbit  have  been  observed  to  charge  to 
thousands  of  volts.  A  joint  program  between  the  Air  Force  and  NASA 
investigated  the  phenomena.  The  program  included  both  theory  and 
experiment.  The  largest  component  was  P78-2,  the  SCATHA  (Spacecraft 
Charging  At  High  Altitudes)  satellite  which  was  instrumented  fully  to  investigate 
magnetospheric  charging.  SCATHA  provided  unrefuteable  evidence  that 
surface  charging  on  spacecraft  can  cause  discharges  and  operational 
anomalies.  As  reported  by  Koons  et.  al.1 ,  on  September  22,  1982  potential 
differences  of  more  than  9500  volts  were  measured  on  the  satellite.  At  the 
same  time,  29  pulses  were  detected  by  the  Transient  Pulse  Monitor.  Seventeen 
of  the  pulses  exceeded  the  maximum  instrument  level  of  7.4  volts.  Coincident 
with  the  discharges  were  three  anomalies,  the  most  severe  resulted  in  a  two 
minute  loss  of  data. 

The  theory  component  of  the  joint  Air  Force  -  NASA  program  culminated  with 
NASCAP,  the  NASA  Charging  Analyzer  Program2,3.  NASCAP  successfully 
modeled  both  frame  and  differential  charging  events  observed  on  SCATHA4.5. 
The  following  briefly  summarizes  the  mechanisms  controlling  spacecraft 
charging  as  they  have  been  included  in  NASCAP.  The  understanding  obtained 
by  developing  NASCAP  lead  to  the  prediction  in  1980  that  large  satellites  in 
polar  orbit  would  charge  from  several  hundred  to  a  few  thousand  volts6.7 .  a 
computer  model,  POLAR  (Potentials  Of  Large  spacecraft  in  Auroral  Region)8, 
was  developed  to  investigate  the  physics  of  auroral  charging  more  fully.  The 
differences  between  auroral  and  geosynchronous  charging  are  discussed 
below.  During  1983,  instruments  on  board  the  Defense  Meterological  Support 
Program  satellite  7  (DMSP  7)  observed  charging  up  to  800  volts9.  Calculations 
using  POLAR  are  compared  with  observed  DMSP  charging.  The  observed 
charging  is  in  agreement  with  the  mechanisms  theorized  in  1980. 

Theory 

The  most  basic  instrument  used  in  laboratory  plasma  experiments  is  the 
Langmuir  probe.  Typically  it  is  a  small  metal  sphere  whose  potential  is  swept 
through  a  limited  range  of  voltages  and  the  current  to  the  sphere  is  measured. 
The  current  is  due  to  charged  particles  from  the  plasma  impinging  upon  the 
sphere.  When  the  sphere  potential  is  large  compared  with  the  kinetic  energy  of 
the  plasma,  only  electrons  are  collected.  When  the  sphere  potential  is  very 
negative,  only  ions  are  collected.  Between  these  two  extremes,  there  is  a 
potential  at  which  the  ion  current  exactly  balances  the  electron  current,  so  that 
the  current  to  the  sphere  is  exactly  zero.  This  potential,  at  which  the  net  current 
is  zero,  is  called  the  floating  potential.  Because  at  a  given  temperature 
electrons  move  so  rapidly  compared  with  ions,  the  floating  potential  is  normally 

negative  a  few  6e,  the  plasma  electron  temperature.  If  the  wire  to  the  probe  is 
cut,  the  probe  rapidly  achieves  the  floating  potential. 
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The  physics  of  spacecraft  charging  can  be  understood  by  regarding  the 
spacecraft  as  a  probe  in  its  local,  ionospheric  plasma.  Since  there  is  no  way 
for  a  continuous  current  to  flow,  the  plasma  particles  rapidly  charge  the 
spacecraft  to  a  few  0e.  The  difference  between  the  laboratory  Langmuir  probe 
and  a  spacecraft  immersed  in  a  magnetospheric  substorm  is  that  the  electron 
energies  are  a  few  volts  in  the  laboratory  and  can  be  tens  of  thousands  of  volts 
in  space.  Laboratory  floating  potentials  are  typically  negative  a  few  volts;  in 
space  potentials  as  high  as  -19,000  volts  have  been  observed10. 


photo  electrons 


Figure  1.  Highly  negative  potentials  can  result  from  the  accumulation  of 
charge  on  spacecraft  surfaces. 


Spacecraft  are  designed  for  purposes  other  than  acting  as  plasma  probes. 
Consequently,  the  interpretation  and  prediction  of  the  spacecraft  potential  are 
complicated  due  to  the  complicated  geometry,  multiple  surface  materials,  and 
the  absence  of  an  easily  accessible  reference  ground.  Each  insulating 
spacecraft  surface  interacts  separately  with  the  plasma  and  is  capacitively  and 
resistively  coupled  to  the  frame  and  other  surfaces.  Rather  than  a  single 
floating  potential,  there  can  be  a  different  one  associated  with  each  surface. 
Computing  surface  potentials  for  a  spacecraft  is  a  considerably  more  complex 
problem  than  computing  the  potential  on  a  conducting,  spherical  probe. 

Not  only  is  additional  complexity  introduced  due  to  geometry  and  insulating 
surfaces,  but  currents  other  than  incident  electrons  and  ions  must  be  included. 
Kilovolt  electrons  generate  secondary  electrons  and  can  be  backscattered  from 
surfacestO.1 1 .  Kilovolt  ions  can  also  generate  secondary  electrons.  The 
current  density  of  low  energy  electrons  generated  by  solar  UV  emission  is  much 
greater  than  the  natural  charging  currents.  Because  of  this,  most  spacecraft 
charging  has  been  observed  during  eclipse,  when  the  spacecraft  is  in  the 
shadow  of  the  earth. 

For  magnetospheric  charging,  the  spacecraft  is  small  compared  with  the  plasma 
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Debye  length.  The  electrons  associated  with  charging  typically  penetrate  less 
than  a  micron  into  the  spacecraft  skin.  Because  of  this,  surface  coatings  play  a 
large  role  in  determining  spacecraft  potentials.  While  the  time  to  achieve 
overall  current  balance  in  very  short,  the  order  of  a  millisecond,  the  time  for 
each  surface  to  achieve  its  own  equilibrium  potential  is  thousands  of  times 
longer.  The  differences  between  the  equilibrium  potentials  of  different  surfaces 
is  referred  to  as  differential  charging. 

n j  =  ne  =  1 06  m-3 
0(  =  0e  ~  1-  30  KeV 
XD  =  1 000  m 

fp  -  105  sec1 

je  *  10"5  amps  m-2 

iPhoto”5><10~5amPs  m_2 

electron  stopping  distance  =  1-10x1 0“7m 

^uniform  charging  sec 

^differential  charging  “  1  ~  1  00  SeC 

Table  1 .  Parameters  typical  of  spacecraft  charging  in  the  magnetosphere. 


In  the  plasmasphere,  spacecraft  usually  float  within  a  few  volts  of  plasma 
potential.  Upon  encountering  a  magnetospheric  substorm,  the  incident  electron 
current  exceeds  the  ion  current  and  the  vehicle  charges  negatively.  As  the 
potential  becomes  negative,  the  electron  current  diminishes  because  not  all  the 
electrons  have  the  energy  to  overcome  the  potential.  If  the  plasma  were  in 
thermal  equilibrium  and  had  a  Maxwellian  distribution  of  energies,  the  electron 
current  would  decrease  exponentially  with  the  negative  potential.  The  spectrum 
of  the  electrons  hitting  the  spacecraft  would  remain  unchanged  since  the  tail  of 
a  Maxwellian  is  still  a  Maxwellian. 

Additional  ions  are  attracted  to  the  spacecraft  as  the  potential  becomes  more 
negative.  For  the  very  low  density  plasma  in  the  magnetosphere,  angular 
momentum  limits  the  collection  of  ions.  The  maximum  impact  parameter  from 
which  ions  are  collected  is  that  for  which  the  ion's  collected  velocity  must  be 
tangent  to  the  spacecraft  in  order  to  conserve  energy  (Figure  2). 

The  balancing  of  ion  and  electron  currents  predicts  a  floating  potential  on  the 
order  of  a  few  times  the  plasma  temperature.  Since  the  electron  current 
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diminishes  exponentially  and  the  ion  current  increases  linearly,  the  principle 
effect  of  the  potential  is  to  decrease  the  electron  current. 


FOR  A  MAXWELLIAN  PLASMA 
n  _  =  density 

0  =  temperature 

REPELLED  SPECIES  IS  ENERGY  LIMITED 

i  =  i0  6  9  for  electrons 

ATTRACTED  SPECIES  IS  ANGULAR 
MOMENTUM  LIMITED 

|  =  |  (A  —  —  \  for  Ions 
1  1  ion' 1  0  * 

POTENTIAL  FOR  CURRENT  BALANCE 

f\ ,  >  .  /  m  > 

CM0  ln(-^y£-iein(^ 


Figure  2.  Orbit  limited  collection  is  based  on  conservation  of  angular 
momentum. 

The  interactions  of  the  incident  electrons  and  ions  with  the  spacecraft  surfaces 
have  a  profound  effect  on  floating  potentials.  The  most  important  process  is 
secondary  electron  emission  1 1 .  Because  secondary  electron  yields  are  so  high 
for  many  surface  materials,  the  spacecraft  floating  potential  is  often  positive!  For 
electrons  with  energies  between  50  eV  and  a  few  Kilovolts,  more  than  one 
secondary  electron  is  emitted  for  every  incident  electron.  This  results  in  a 
positive  charging  current.  Only  when  the  electron  energies  exceed  several 
thousand  volts  will  the  spacecraft  charge  negatively.  Backscatter  yields  are  less 
than  unity  and  vary  little  with  energy.  Ion  generated  secondary  electrons 
enhance  the  ion  current  and  act  to  reduce  overall  charging  levels. 

A  kapton  sphere  immersed  in  a  5  keV,  106  rn-3  plasma  charges  to  -3640  volts. 
The  electron  current  which  drives  the  charging  is  less  than  2X10-6  A/m2,  more 
than  half  of  which  is  immediately  cancelled  by  secondaries  and  backscatter. 

Ion  generated  secondary  electrons  effectively  triple  the  incident  ion  current. 
Equilibrium  occurs  at  a  potential  less  than  the  plasma  temperature.  Because  of 
the  secondaries  and  backscatter,  current  balance  is  effected  equally  by 
diminishing  the  electron  and  increasing  the  ions  currents.  Because  the  incident 
electron  spectrum  remains  Maxwellian,  electron  generated  secondaries  and 
backscattered  electrons  remain  a  constant  fraction  of  the  incident  current  as 
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the  spacecraft  charges.  The  ion  generated  secondaries  increase  compared 
with  the  incident  ion  current  because  the  energy  of  the  ions  increases  as  the 
spacecraft  potential  becomes  more  negative.  Ion  generated  secondary 
electron  yields  peak  for  ion  energies  of  several  tens  of  kilovolts. 


energy  (eV) 

Figure  3.  Electron  generated  secondary  electron,  backscatter  and  proton 
generated  secondary  electron  yields  for  kapton. 

Photoelectron  current  densities  are  about  10  4  A/m2,  typically  an  order  of 
magnitude  greater  than  incident  electron  currents,  even  before  secondaries  are 
taken  into  account.  As  a  result,  in  sunlight  high  negative  potentials  are  rarely 
observed  on  spacecraft.  Charging  in  sunlight  occurs  when  the  low  energy 
photoelectrons  can't  escape  from  the  spacecraft  because  of  potential 
barriers1 2. 13 .  in  an  intense  substorm,  spacecraft  surfaces  shadowed  from  the 
sun  slowly  charge  to  thousands  of  volts  negative,  while  the  sunlit  surfaces 
remain  a  few  volts  positive.  After  some  time,  a  saddle  point  in  the  potential 
develops  in  front  of  the  sunlit  surfaces.  Because  the  saddle  point  is  driven  by 
the  photoelectron  kinetic  energy,  it’s  height  can  never  exceed  a  few  volts.  All 
the  surfaces  then  charge  negative  at  a  rate  corresponding  with  differential 
charging,  typically  a  few  hundred  volts  per  minute.  Sunlight  charging  is  a  multi 
dimensional  effect  calculable  only  by  models,  such  as  NASCAP  and  POLAR, 
which  solve  Poisson’s  equation  in  three  dimensions14. 

Auroral  Charging 

The  model  for  enhanced  collection  described  above  is  generally  referred  to  as 
orbit  limited  current  collection.  It  is  appropriate  when  the  potential  has  a  range 
larger  than  the  largest  impact  parameter  and  is  sufficiently  well  behaved  so  that 
no  angular  momentum  barriers  exist.  Potentials  that  vary  mere  slowly  than  with 
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the  inverse  of  the  radius  squared  satisfy  these  conditions.  At  geosynchronous 
orbit,  the  plasma  is  so  thin  that  little  shielding  occurs  and  the  spacecraft 
potential  drops  roughly  as  the  inverse  of  the  radius. 

Orbit  limited  collection  is  not  valid  at  lower  altitudes,  where  the  plasma  is  much 
denser.  If  the  spacecraft  were  to  charge  negative,  the  additional  ions  collected 
would  shield  and  thus  limit  the  range  of  the  potential.  When  the  collection 
radius  is  limited  by  shielding  due  to  the  charge  of  the  collected  species,  space 
charge  limited  collection  is  the  appropriate  theory.  To  model  spacecraft 
charging  in  ionospheric  plasma  with  densities  greater  than  about  109  m-3, 
space  charge  limited  collection  models  must  be  used. 


Figure  4.  Current  vs  voltage  for  a  kapton  sphere  in  a  5  keV,  106  m-3 
Maxwellian  plasma.  The  floating  potential  is  -3640  volts. 


The  energetic  electrons  which  charge  spacecraft  in  low  polar  orbit  are  the  same 
that  generate  the  aurora  borealis.  While  of  similar  origin  to  substorm  electrons 
in  the  magnetosphere,  the  auroral  electron  fluxes  can  be  as  much  as  a  hundred 
times  as  intense.  Some  of  the  enhanced  intensity  comes  from  the  convergence 
of  the  magnetic  field  lines  as  they  approach  the  poles.  Measured  fluxes  can  be 
fit  well  using  the  analytical  form  suggested  by  Fontheim  et.  al.15. 

Charging  of  large  objects  in  polar  orbit  is  determined  by  the  balance  of  the  net 
auroral  flux  and  the  space  charge  limited  ion  flux.  That  much  higher  spacecraft 
potentials  are  obtained  due  to  space  charge  limiting  was  first  predicted  in 
19806.  A  comparison  between  the  potentials  obtained  ignoring  space  charge 
with  that  including  it  shows  more  than  an  order  of  magnitude  difference.  The 
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space  charge  limited  result  agrees  with  observation. 

Spacecraft  Charging  Computer  Codes 

The  NASA  Charging  Analyzer  Erogram  (NASCAP),  was  the  first  three 
dimensional  spacecraft-environment  interaction  code.  NASCAP  solves 
Poisson's  equation  for  spacecraft  potentials  due  to  the  accumulation  of  charge 
from  the  magnetospheric  plasma.  Asymmetries  are  introduced  by  potential 
differences  between  various  spacecraft  surfaces  due  to  differences  in  material 
response  and  solar  illumination.  The  surface  potential  variations  create  saddle 
point  barriers  which  control  low  energy  secondary  currents.  NASCAP 
represents  the  spacecraft  using  a  restricted  constructive  geometry  technique 
with  building  blocks  of  cubes  and  segments  of  cubes.  Because  the  potentials 
vary  on  a  short  distance  scale  near  the  spacecraft  and  more  slowly  far  away, 
NASCAP  zones  the  space  exterior  to  the  satellite  using  a  sequence  of  nested 
finite  element  grids.  The  NASCAP  code  was  first  used  to  model  the  SCATHA 
satellite.  Since  then,  NASCAP  has  been  used  throughout  the  world  for 
spacecraft  charging  analysis  of  geosynchronous  satellites. 

The  Air  Force  Geophysics  Laboratory  developed  POLAR  code  models  the 
charging  of  large,  polar  orbiting  spacecraft.  POLAR  also  calculates  the  plasma 
wake  as  part  of  the  charging  analysis.  POLAR  zones  space  using  a  staggered 
grid  so  that  most  of  the  computational  space  is  aligned  along  the  flow  direction. 
The  steady  state  sheath  structure  is  determined  by  iterating  between  Poisson’s 
equation  for  potentials  and  particle  tracking  for  space  charge  density.  The 
POLAR  wake  model  has  been  compared  with  in  situ  measurements  of  the  wake 
created  by  the  Space  Shuttle  Orbiterl 6. 

Calculations  of  auroral  electron  induced  charging  of  one  of  the  Defense 
Meterological  Support  Program  (DMSP)  satellites  have  recently  been 
compared  with  measurement.  With  the  assistance  of  M.  S.  Gussenhoven,  one 
of  the  observed  charging  electron  spectra  from  December  1 983  was  fit  to  the 
Fontheim  form  for  use  in  POLAR.  The  POLAR  calculation  resulted  in  a 
spacecraft  potential  of  -230  volts,  in  close  agreement  with  the  observed 
potential  of  -215  volts.  The  differences  are  smaller  than  the  possible  sources  of 
error  in  such  a  calculation.  The  incident  electron  spectrum  is  measured  using 
just  a  few  bins,  the  spacecraft  surface  properties  change  with  time  in  orbit,  the 
details  of  the  ionic  composition  are  not  known  that  accurately.  In  general, 
differences  larger  than  the  10%  shown  in  this  particular  case  are  to  be  expected 
between  calculations  and  observations. 

Conclusion 

Spacecraft  charging  is  important  because  it  has  been  shown  to  cause 
operational  upsets  on  satellites.  While  of  limited  importance  for  small  satellites 
in  low  polar  orbit,  recent  data  confirms  predictions  that  charging  will  become 
more  important  as  the  size  of  polar  orbiting  spacecraft  increase. 
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ABSTRACT 

The  three-dimensional  computer  codes  NASCAP/LEO  and  POLAR  have  been  used 
to  calculate  current  collection  by  the  mother  payload  of  the  Charge-2  rocket  under  condi¬ 
tions  of  positive  and  negative  potential  up  to  several  hundred  volts.  For  negative  bias 
(ion  collection)  the  calculations  fall  about  twenty-five  percent  above  the  data,  indicating 
that  the  ions  were  less  dense,  colder,  or  heavier  than  the  input  parameters.  For  positive 
bias  (electron  collection)  NASCAP/LEO  and  POLAR  calculations  show  similar  agree¬ 
ment  with  the  measurements  at  the  highest  altitudes.  This  indicates  that  the  current  is 
classically  magnetically  limited,  even  during  electron  beam  emission.  However,  the  cal¬ 
culated  values  fall  well  below  the  data  at  lower  altitudes.  We  suggest  that  beam-plasma- 
neutra!  interactions  are  responsible  for  the  high  values  of  collected  current  at  altitudes 
below  240  kilometers. 
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I.  Introduction 

For  many  years  electron  beams  have  been  operated  from  spacecraft  and  rockets  to 
study  long  distance  transport  of  electrons  along  the  earth’s  magnetic  field  (Winckler, 
1980).  These  studies  were  often  compromised  by  charging  of  the  spacecraft  and/or 
severe  perturbation  of  its  environment  as  a  result  of  beam  operations.  The  study  of  these 
interactions  has  proved  to  be  a  complex  and  interesting  field  in  itself. 

To  analyze  a  beam  charging  event,  one  first  asks  whether  the  observations  can  be 
explained  by  simple,  first-order,  passive  current  collection.  This  turns  out  to  be  a 
difficult  question,  rnce  for  the  relevant  conditions  (intermediate  thickness  sheath  with 
magnetic  field)  analytic  theories  exist  only  for  spherical  spacecraft.  In  this  paper  we 
describe  the  application  of  the  three-dimensional  computer  codes  NASCAP/LEO  and 
POLAR  to  calculate  charge  collection  by  the  highly  cylindrical  mother  vehicle  of  the 
Charge-2  rocket.  Calcu.arions  and  measurements  for  negative  bias  (ion  coiiection)  have 
been  reported  elsewhere  (Neubert  et  al.,  1989),  and  will  only  be  briefly  summarized  here. 
The  measurements  for  positive  bias  have  also  been  reported  previously  (Myers  et  al., 
1989).  In  that  paper  the  difficulty  of  assessing  the  results  using  theories  which  apply 
only  to  spheres  and  infinite  cylinders  is  made  apparent.  The  contribution  of  this  paper  is 
the  interpretation  of  the  positive  bias  results  through  the  use  of  codes  which  correctly 
account  for  the  true  three-dimensional  geometry  of  this  problem. 

1.1.  The  Charge-2  Rocket 

The  Charge-2  payload  was  launched  on  a  Black  Brant  VB  from  White  Sands  Mis¬ 
sile  Range  in  New  Mexico  in  December  1985.  The  payload  consisted  of  mother  and 
daughter  sections  which  were  electrically  connected  by  an  insulated  tether.  The  mother 
payload  was  a  cylinder  of  diameter  0.44  m  and  length  3.26  m;  the  0.41  m  nose  section 
was  coned  down  to  a  diameter  of  0.33  m;  the  total  exposed  area  was  4.68  m2.  The 
daughter  payload  was  a  cylinder  of  diameter  0.44  m  and  length  1.62  m. 

Negative  mother  potentials  were  achieved  by  applying  a  potential  between  the 
mother  chassis  and  the  tether  end.  The  daughter  potential  required  to  collect  enough 
electrons  to  balance  the  ion  current  to  the  mother  is  negligible,  so  that  the  mother  poten¬ 
tial  is  equal  to  the  applied  potential.  To  maintain  a  steady  state,  the  ion  current  collected 
by  the  mother  vehicle  must  be  equal  to  the  tether  current,  which  is  easily  measured. 

Positive  mother  potentials  were  achieved  by  means  of  electron  beam  operations. 
The  electron  current  collected  by  the  mother  is  then  equal  to  the  beam  current  less  the 
tether  current. 

The  mother  payload  also  contained  an  array  of  four  "PLP"  probes  located  on  a 
boom  at  locations  25,  50,  75,  and  100  cm  from  the  rocket  surface.  The  intent  was  to  use 

g 

these  to  measure  the  sheath  profile.  The  10  ft  probe  impedance  was  far  too  low  to  do 
this  in  the  negative  bias  case,  but  should  have  worked  well  for  positive  bias  where  the 
current  levels  are  higher.  The  rocket  potential  under  positive  bias  conditions  was 
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inferred  from  the  PLP  probes. 

1.2.  The  NASCAP/LEO  Computer  Code 

NASCAP/LEO  (Mandell  et  al„  1982)  is  a  computer  code  developed  by  the  S- 
CUBED  Division  of  Maxwell  Laboratories  under  contract  to  NASA/Lewis  Research 
Center  to  model  the  electrostatic  interactions  of  a  high  voltage  spacecraft  with  a  plasma 
under  conditions  commonly  encountered  in  low  earth  orbit.  Because  NASCAP/LEO  is 
fully  three-dimensional,  it  is  able  not  only  to  treat  the  actual  spacecraft  shape,  but  also  to 
account  for  asymmetries  introduced  by  spacecraft  motion  and  magnetic  field  effects. 

NASCAP/LEO  accepts  as  a  geometrical  specification  a  finite  element  surface  model 
produced  by  any  of  several  standard  CAD  programs.  This  model  is  then  embedded  in  a 
cubic  grid.  Within  this  "primary  grid",  locally  enhanced  resolution  near  small  but  impor¬ 
tant  spacecraft  features  may  be  achieved  through  the  use  of  cell  subdivision.  To  include 
a  large  volume  of  space  about  the  spacecraft,  the  primary  grid  may  be  nested  within  suc¬ 
cessive  "outer  grids".  Each  outer  grid  has  the  same  number  of  grid  points  as  the  next 
inner  grid,  but  rwice  the  physical  spacing  between  the  grid  points,  so  that  it  includes  eight 
times  the  volume. 

NASCAP/LEO  calculates  electrostatic  potentials  in  the  plasma  surrounding  a  space¬ 
craft  by  using  finite  element  and  conjugate  gradient  methods  to  solve  the  variational  form 
of  Poisson’s  equation.  The  space  charge  appearing  in  Poisson’s  equation  is  represented 
as  an  analytic  function  of  the  plasma  density  and  temperature  and  the  local  potential  and 
electric  field.  This  function  reduces  to  linear  screening  at  low  potentials,  and  takes 
account  of  particle  acceleration  and  convergence  effects  at  high  potentials. 
(NASCAP/LEO  does  not  account  for  modification  of  space  charge  and  potential  structure 
due  to  magnetic  fields.)  Boundary  conditions  at  object  surfaces  may  be  either  specified 
potential  or  specified  electric  field. 

Currents  to  spacecraft  are  calculated  using  the  "sharp  sheath  edge"  approximation. 
A  specified  electrostatic  potential  contour  is  designated  as  the  sheath  surface.  Macropar¬ 
ticles  are  generated  to  represent  the  plasma  thermal  current  (modified  by  spacecraft 
motion  effects)  passing  through  elements  of  area  of  the  sheath  surface.  Each  macroparti¬ 
cle  is  then  tracked  (in  the  calculated  electric  fields  and  specified  magnetic  field)  until  it 
impinges  on  the  spacecraft  or  some  other  termination  condition  (e.g.,  escaping  the  grid) 
occurs.  Thus  the  total  current  to  the  spacecraft  and  the  current  distribution  over  the 
spacecraft  surface  is  determined. 

NASCAP/LEO  also  contains  surface  charging  algorithms,  solar  array  models,  and 
hydrodynamic  ion  models  which  are  not  relevant  to  the  subject  of  this  paper. 
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1.3.  The  POLAR  Computer  Code 

POLAR  (Cooke  ex  al.,  1985)  is  a  computer  code  developed  by  the  S-CUBED  Divi¬ 
sion  of  Maxwell  Laboratories  under  contract  to  Air  Force  Geophysics  Laboratory  to 
model  charging  of  a  large  spacecraft  in  the  polar  auroral  environment.  Thus  POLAR  has 
an  environment  model  capable  of  representing  the  energetic  and  anisotropic  panicle 
fluxes  observed  in  polar  regions.  However,  the  POLAR  calculation  strategy  for  this  work 
did  not  take  advantage  of  the  sophisticated  environment  model.  Rather,  POLAR  was 
used  to  calculate  the  current  collected  by  the  spacecraft  in  a  fully  self-consistent  space 
charge  field,  to  augment  the  NASCAP/LEO  results  which  ignore  magnetic  field  effects 
on  space  charge. 

POLAR  requires  that  spacecraft  be  modeled  in  a  restricted  geometry  set  similar  to 
that  used  in  the  NASCAP/GEO  code.  This  allows  fully  three-dimensional  modeling  for 
spacecraft  models  of  modest  complexity. 

The  POLAR  grid  structure  was  designed  to  encompass  the  wake  region  behind  a 
spacecraft.  Thus,  the  cubic  grid  is  made  up  of  a  sequence  of  "bread  slices"  aligned 
approximately  normal  to  the  spacecraft  velocity.  Because  the  spacecraft  velocity  may 
not  coincide  with  a  coordinate  direction,  the  "bread  slices"  may  be  staggered  in  order  to 
be  centered  about  the  spacecraft  wake.  With  this  structure,  POLAR  has  neither  subdi¬ 
vided  grids  nor  outer  grids,  but  maintains  constant  resolution  throughout  the  computa¬ 
tional  space.  Since  there  is  no  loss  of  resolution  at  the  sheath  surface,  the  calculated 
current  to  a  spacecraft  has  less  dependence  on  the  choice  of  sheath  boundary  potential 
than  is  the  case  with  NASCAP/LEO. 

The  most  important  POLAR  capability  for  the  present  purpose  is  its  ability  to  use 
particle  tracking  results  to  modify  the  space  charge  calculated  by  analytic  approxima¬ 
tions.  Thus,  starting  from  a  potential  solution  based  on  analytic  space  charge  as  in 
NASCAP/LEO,  POLAR  can  take  account  of  the  modification  of  particle  trajectories  by 
magnetic  fields  to  calculate  self-consistent  space  charge  and  potential  distributions. 

2.  Calculations  for  Negative  Bias 

Comparison  of  NASCAP/LEO  results  to  the  measured  current  collected  by  the 
Charge-2  mother  payload  under  negative  bias  conditions  has  been  published  elsewhere 
(Neubert  et  al.,  1989).  In  addition  to  calculating  the  total  current  to  the  rocket, 
NASCAP/LEO  did  a  fairly  good  job  on  the  more  difficult  problem  of  calculating  the  bias 
on  the  PLP  probes  operating  in  their  "floating"  mode.  Figure  1  shows  the  NASCAP/LEO 
model  of  Charge-2  in  the  grid,  with  enhanced  resolution  in  the  neighborhood  of  the  PLP 
booms. 

Here  we  will  take  a  somewhat  different  point  of  view.  By  forcing  NASCAP/LEO  to 
calculate  the  correct  current  for  an  equivalent  spherical  model  (Langmuir  and  Blodgett, 
1924),  we  can  have  confidence  that  the  difference  in  calculated  current  between  the 
sphere  and  Charge-2  is  due  to  the  effects  of  non-spherical  geometry  and  spacecraft 
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motion  and  not  from  poor  resolution  at  the  sheath  boundary.  The  calculated  results  will 
be  compared  with  the  measurements  as  well.  All  calculations  in  this  section  are  done 
with  the  plasma  parameters  and  computational  parameters  shown  in  Table  1. 


Table  1:  Plasma  and  Computational  Parameters 


Plasma  Density 

8xl010m'3 

Plasma  Temperature 

0.05  eV 

Debye  Length 

0.588  cm. 

Electron  Thermal  Current 

4.8x1  O'4  A-m'2 

Ion  [O^]  Thermal  Current 

2. 8x10^  A-m2 

Magnetic  Field 

0 

Spacecraft  Velocity 

0 ;  580  m-sec'1 

Primary  Grid  Resolution 

10.8  cm. 

Sheath  Boundary  Potential 

(varies) 

Charge-2  Model  Area 

4.63  m2 

Sphere  Model  Area 

4.45  m2 

The  major  source  of  uncertainty  in  NASCAP/LEO  current  calculations  lies  in  the 
choice  of  sheath  boundary  potential.  In  the  common  case  that  the  plasma  temperature  is 
very  low  and  many  Debye  lengths  are  contained  in  a  computational  zone,  the  optimal 
sheath  boundary  potential  will  be  well  above  the  plasma  temperature.  The  reason  for  this 
follows  from  the  theory  of  space  charge  limited  current  flow.  For  a  plasma  of  density  n 
(m  3]  and  temperature  0  [eV]  the  potential  a  distance  Ax  in  from  the  physical  sheath  will 
be 

<|>(Ax)  =  5-lxlfr6  [(Ax)4  0  n2],/3 . 

For  the  parameters  in  Table  1,  this  value  is  1.8  volts  in  the  primary  grid,  4.5  volts  in  the 
first  outer  grid,  and  1 1.4  volts  in  the  second  outer  grid.  However,  NASCAP/LEO  allows 
at  most  an  order-of  magnitude  potential  drop  per  zone  due  to  space  charge  effects.  Thus, 
a  minimum  appropriate  choice  for  the  sheath  boundary  potential  is  one  order-of- 
magnitude  below  0(Ax),  which  is  certainly  well  above  the  plasma  temperature  of  .05 
volts.  Experience  has  shown  that  a  poor  choice  of  sheath  boundary  potential  can  lead  to 
errors  of  as  large  as  thirty  percent  in  the  calculation  of  collected  current. 

In  this  work  we  have  taken  care  to  minimize  the  uncertainty  in  the  choice  of  sheath 
boundary  potential.  We  have  determined  the  sheath  boundary  potential  required  to 
obtain  the  correct  (Langmuir  and  Blodgett,  1924)  current  to  a  stationary  sphere  of  area 
similar  to  the  Charge-2  model  with  the  same  plasma  and  computational  conditions,  and 
used  that  sheath  boundary  potential  to  calculate  the  current  to  the  Charge-2  model.  The 
sheath  boundary  potential  values  used  ranged  from  1.0  volts  for  the  -10  volt  cases  to  3.4 
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volts  for  the  -400  volt  cases.  By  this  procedure,  we  believe  the  computational  error  asso¬ 
ciated  with  choice  of  sheath  boundary  potential  has  been  reduced  to  a  few  percent  at 
most. 

Table  2  shows  the  exact  current  to  a  sphere  of  area  4.45  m  and  the  calculated 
currents  to  the  Charge-2  model,  using  the  same  sheath  boundary  potential  and  taking 
account  of  the  spacecraft  velocity  of  580  m-sec‘l.  The  effect  of  spacecraft  motion  is  to 
increase  the  current  by  -  40  percent  at  low  potential,  dropping  to  -  15  percent  at  high 
potential.  We  have  also  enhanced  the  collected  current  at  -400  volts  by  fifty  percent  of 
the  calculated  incident  ion  current  to  account  for  secondary  electron  generation  by  0+ 
impact.  This  process  has  a  threshold  at  about  200  eV,  and  a  yield  of  about  one-half  at 
400  eV  (Dietz  and  Sheffield,  1975). 


Table  2.  Calculated  Currents  [pA] 
to  NASCAP/LEO  Sphere  and  Charge-2  Models 


Potential 

Sphere  Model 
v  =  0 

Charge-2  Model 
v  =  580  m  sec"1 

-10 

28.9 

52 

-20 

38.8 

68 

-50 

62.2 

98 

-100 

94 

140 

-200 

148 

203 

-400 

239 

457a 

a.  Includes  fifty  percent  enhancement  due  to  secondary  electron  emission. 

On  figure  2  the  calculations  fall  above  the  data  by  about  twenty-five  percent.  Part 
of  the  discrepancy  is  due  to  the  presence  of  a  substantial  NO+  component  in  the  ion 
population  (Wahl,  1988;  Rawer  and  Bradley,  1987),  which  lowers  the  ion  thermal 
current;  a  25  percent  NO+  component  would  reduce  the  collected  current  by  about  ten 
percent.  The  remaining  discrepancy  must  be  attributed  to  errors  in  the  ion  density  value 
used  in  the  calculation  or  to  other  physical  uncertainties.  Note  that  Neubert  et  al.  (1989) 
used  an  ion  density  of  4x10  m  and  a  temperature  of  0.1  eV  for  all  calculations. 

3.  Calculations  for  Positive  Bias 

In  the  case  of  positive  bias  (electron  collection)  the  earth’s  magnetic  field  plays  a 
dominant  role  in  limiting  the  current  collected  by  a  spacecraft.  Parker  and  Murphy 
(1967)  showed,  using  conservation  of  energy  and  angular  momentum,  that  a  bound  on 
the  current  collected  by  a  sphere  (actually,  by  a  surface  of  revolution  about  an  axis  paral¬ 
lel  to  the  magnetic  field)  is  given  by 

I  <  2na2  [  1  +  (8e<t>/mo)2a2),/4  ] 
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where  a  is  the  sphere  radius  and  to  is  the  electron  gyrofrequency,  eB/m.  As  the  Charge-2 
cylinder  axis  was  normal  to  the  magnetic  field,  angular  momentum  is  not  conserved  so 
that  the  Parker-Murphy  bound  does  not  strictly  apply.  The  bound  will  also  not  apply  in 
the  presence  of  ionization,  turbulence,  or  other  collisional  effects.  Nonetheless,  it  is  a 
useful  guide  to  whether  the  earth’s  magnetic  field  is  the  dominant  current  limiting 
mechanism. 

Table  3a  gives  the  environmental  parameters  and  spacecraft  potential  for  the  previ¬ 
ously  published  (Myers  et  al.,  1989)  current  measurements,  and  Table  3b  gives  the  meas¬ 
ured  and  calculated  currents.  In  all  cases,  the  magnetic  field  was  taken  to  be  0.4  gauss, 
the  plasma  temperature  0.052  eV,  and  the  NASCAP/LEO  sheath  boundary  potential  3.8 
volts  (chosen  to  give  the  correct  sheath-limited  current  to  the  equivalent  sphere  in  the 
260  km  environment).  The  current  collected  by  the  spacecraft  is  magnetically  limited,  as 
most  of  the  electrons  entering  the  sheath  cannot  cross  enough  field  lines  to  be  collected. 

For  the  four  data  points  above  240  kilometers,  the  measured  and  calculated  currents 
are  in  excellent  agreement.  Note  that  both  exceed  the  Parker-Murphy  bound,  indicating 
the  difficulty  of  applying  a  symmetry-requiring  theory  to  a  non-symmetric  object.  (We 
would  expect  a  similar,  though  smaller,  effect  in  calculating  magnetically  limited  current 
to  the  daughter  rocket,  so  that  the  daughter  current  under  positive  bias,  shown  in  Myers  et 
al.  (1989)  as  lying  on  the  Parker-Murphy  bound  for  the  equivalent  sphere,  is  actually 
slightly  below  the  correct  magnetically  limited  current.)  These  data  confirm  the  role  of 
the  magnetic  field  in  limiting  the  collected  current,  as  the  measurements  fall  far  below 
the  sheath-limited  current  (i.e.,  current  limited  only  by  space  charge  effects).  The  calcu¬ 
lation  of  magnetically  limited  current  is  fairly  insensitive  to  the  choice  of  sheath  boun¬ 
dary  potential;  it  increases  slightly  with  the  sheath  boundary  potential. 
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Table  3a.  Environments  for  Electron  Current  Data 


1  Altitude 

Potential 

Neutral  Density* 

Plasma  Density 

Ikmj 

[volts] 

[10l5m'3] 

[  1010  m'3] 

165 

390 

36 

0.46 

168 

150 

32 

0.48 

232 

475 

6 

5.6 

251a 

560 

3.9 

8.0 

251b 

440 

3.9 

8.5 

256 

440 

3.5 

9.0 

260 

440 

3.3 

9.3 

Table  3b.  Calculated  and  Measured  Collected  Current  [mA] 

Altitude 

Parker- 

Measured 

NASCAP/LEO 

NASCAP/LEO 

POLAR 

[km]  Murphyb 

Collected  Current 

Sheath  Current  Self-Consistent 

165 

0.4 

35.8 

0.6 

8.3 

- 

168 

0.3 

6 

0.5 

3.5 

- 

232 

5.5 

20.4 

7.9 

42 

- 

251a 

8.5 

12.2 

13.0 

60 

- 

251b 

8.1 

14 

12.9 

52 

- 

256 

8.6 

15.6 

13.5 

54 

- 

260 

8.9 

18 

14.1 

55 

22 

a.  MSIS-86. 

b.  Calculated  for  sphere  of  radius  0.6  meters. 

A  POLAR  code  calculation  was  done  for  the  highest  altitude  measurement. 
POLAR,  after  performing  a  calculation  similar  to  NASCAP/LEO,  used  particle  trajectory 
information  to  obtain  space  charge  densities  and  electrostatic  potentials  consistent  with 
the  trajectories  in  the  presence  of  magnetic  field.  The  POLAR  model  of  Charge-2  was  a 
right  octagonal  cylinder  with  length  3.3  meters  and  area  4.85  m  ,  and  the  code  resolution 
was  0.14333  meters.  The  higher,  self-consistent  space  charge  densities  brought  the 
sheath  surface  closer  to  the  object  surface,  resulting  in  a  fifty  percent  increase  in  current 
over  the  non-self-consistent  calculation.  We  expect  that  the  percentage  change  will  be 
less  for  lower  ambient  plasma  densities.  This  self-consistency  effect  explains  why  the 
measured  current  is  consistently  above  the  NASCAP/LEO  calculated  current.  It  is  worth 
noting  that  the  POLAR  calculated  current  of  22  mA  does  fall  below  the  Parker-Murphy 
bound  (30  mA)  for  r  sphere  with  diameter  equal  to  the  rocket  length. 

The  three  data  points  below  240  kilometers  clearly  show  measured  currents  exceed¬ 
ing  the  magnetically  limited  value.  At  232  kilometers,  the  collected  current  is  about  dou¬ 
ble  that  calculated,  but  well  under  the  sheath-limited  current;  at  165  and  168  kilometers 
the  measurement  is  one  to  two  orders  of  magnitude  above  the  calculation.  Simple 
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ionization  models  are  unable  to  predict  any  substantial  effect  within  the  sheath.  We  are 
forced  to  conclude  that  a  complex  beam-plasma  interaction  is  taking  place  at  these  lower 
altitudes. 

4.  Conclusions 

Through  the  use  of  three-dimensional  modeling  we  are  able  to  make  unambiguous 
comparison  with  measurements,  in  contrast  to  with  theories  valid  only  for  symmetric 
geometry. 

It  has  been  published  previously  that  NASCAP/LEO  is  able  to  obtain  excellent 
agreement  between  calculated  and  measured  currents  under  conditions  of  negative  bias. 
The  effects  of  non-spherical  geometry  and  of  spacecraft  motion  can  be  separated.  Also, 
NASCAP/LEO  was  able  to  obtain  good  values  for  the  biases  on  the  PLP  probes. 

Because  NASCAP/LEO  has  poor  resolution  at  the  sheath  surface,  an  intelligent 
choice  of  sheath  boundary  potential  is  necessary  to  obtain  accurate  current  collection 
values.  A  non-optimal  sheath  boundary  potential  will  lead  to  errors  of  a  few  tens  of  per¬ 
cent.  By  comparison  with  analytic  models  in  the  same  computational  environment  we 
obtained  computational  results  with  probable  accuracy  of  better  than  ten  percent.  By 
comparing  the  computational  results  for  nominal  ion  parameters  with  the  measurements, 
we  conclude  that  errors  are  due  to  the  presence  of  an  NO+  ion  component  and  probably  a 
lower  ion  density  than  the  reported  value.  The  onset  of  secondary  electron  emission 
above  200  eV  must  be  taken  into  account. 

Under  positively  biased  (by  electron  beam)  conditions,  NASCAP/LEO  calculations 
of  magnetically  limited  current  are  in  good  agreement  with  previously  published  experi¬ 
mental  results  at  altitudes  over  240  kilometers.  As  the  measured  currents  exceed  the 
spherical  Parker-Murphy  bound,  it  was  not  clear  before  these  calculations  whether  simple 
magnetically  limited  collection  is  adequate  to  account  for  the  measured  currents.  These 
calculations  demonstrate  unambiguously  that  magnetically  limited  collection  was 
observed  above  240  kilometers.  Other  mechanisms  are  needed  to  explain  the  observa¬ 
tions  at  lower  altitudes. 

As  NASCAP/LEO  does  not  treat  the  modification  of  space  charge  by  magnetic  field 
effects,  we  have  run  the  POLAR  code  to  study  the  effect  of  fully  self-consistent  space 
charge  on  the  highest  altitude  measurement  For  this  case,  fully  self-consistent  space 
charge  led  to  an  increase  in  current  of  about  fifty  percent,  bringing  the  calculation  into 
better  agreement  with  the  measurement,  and  explaining  the  general  trend  that  the  non¬ 
self-consistent  calculations  lie  below  the  high  altitude  measurements.  As  the  plasma 
density  decreases  and  the  sheath  lies  farther  beyond  the  Parker-Murphy  radius,  we  expect 
the  difference  in  space  charge  treatments  to  have  less  effect  on  the  calculated  currents. 
Note  also  that  the  fully  self-consistent  result  is  about  twenty  percent  high  relative  to  the 
measurement;  this  is  a  comparable  error  to  the  negative  bias  results,  again  suggesting  that 
the  actual  densities  were  below  the  nominal  values. 
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In  summary,  the  calculations  clearly  demonstrate  that  the  published  current  collec¬ 
tion  data  above  240  kilometers  can  be  understood  in  terms  of  simple,  quasi-static  sheath 
prex'es-es  with  magnetic  liuiiting  for  electrons,  t  he  three  low  altitude  electron  current 
measurements  lie  well  above  the  predictions  of  such  theories  and  indicate  that  more  com¬ 
plex  physical  mechanisms  are  at  work. 
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Figure  Captions 

Figure  1.  NASCAn/LEO  gridding  around  the  Charge-2  model,  showing  fine  resolution 
about  the  PLP  booms  and  halving  of  resolution  in  the  outer  grid. 

Figure  2.  Ion  collection  current:  points  and  line:  this  paper;  horizontal  error  bars:  flight 
results. 
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